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Summary 


A detailed exact solution based on the Sanders-Koiter linear shell theory is presented for 
laminated-composite circular cylinders that undergo axisymmetric deformations. The overall 
solution is formulated in a general, systematic way and is obtained from the solution of a single 
fourth-order, nonhomogeneous ordinary differential equation with constant coefficients in which 
the radial displacement is the dependent variable. Moreover, the full response-coupling effects 
associated with generally laminated wall construction are included in the solution. A key element 
of the analysis is the use of the positive-definiteness of the strain energy to define uniquely the 
form of the basis functions spanning the solution space of the ordinary differential equation. 
Loading conditions are considered that include axisymmetric edge loads, surface tractions, and 
temperature fields. Likewise, all possible axisymmetric displacement and traction boundary 
conditions are considered. 

Results are presented for five examples that demonstrate a wide range of behavior for 
specially orthotropic and fully anisotropic cylinders. Several of the results illustrate how 
laminated-wall construction, loading and support conditions, and cylinder length affect the size 
of localized zones of bending deformations that typically occur near the cylinder ends. For some 
of the examples, the results indicate highly localized bending deformations that attenuate rapidly. 
For other examples, the bending deformations exhibit very little attenuation. Approximate 
solutions are also given for several of the examples that add credibility to the corresponding 
results obtained from the exact solution and that provide insight into the possible pitfalls of 
making certain types of simplifications in a preliminary-design process. For example, several 
results are presented that illustrate conditions for which sizing a laminated-composite cylinder 
based on an approximate linear membrane analysis is inappropriate. 

Results obtained from the exact solution presented herein also have the potential to serve as 
fundamental benchmarks for the validation of numerical solution methods such as the finite- 
element method. Likewise, results obtained from the exact solution have the potential to provide 
insight into the extent of discretization that might be needed for analyzing accurately shells with 
zones of localized bending caused by cutouts, stiffeners, or other types of discontinuities and 
gradients. 
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Introduction 

Weight and structural performance are key elements in the design of aircraft and spacecraft. 
For example, weight directly affects fuel costs incurred by airlines, which affects consumers and 
commerce. In addition, structural performance affects the safety of the flying public. Similarly, 
the present cost for placing a spacecraft in Earth orbit ranges from approximately $5,000 to 
$10,000 per pound, depending on the altitude. Thus, structural weight is a significant design 
driver in launch-vehicle structures that operate near the Earth, such as the Space Shuttle external 
tank and solid rocket boosters shown in figure 1 , and in vehicles that operate far from Earth such 
as the conceptual lunar lander shown in figure 2. In recent studies of newer vehicles like a lunar 
lander, fiber-reinforced, laminated-composite materials have received significant attention for use 
in primary structural members. These advanced composite materials have higher stiffness-to- 
density ratios, higher strength-to-density ratios, and very low coefficients of thermal expansion, 
compared to common isotropic materials. Additionally, because laminated composites are 
composed of orthotropic layers, oriented at various angles, the potential exists to tailor a laminate 
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to meet a specific set of high-performance, or special purpose, structural requirements. Thus, 
laminated composites offer a design space that is much wider than that for isotropic materials. 
However, along with the potential structural performance benefits of laminate tailoring comes a 
much more complicated analysis and design process. 

Efforts to reduce structural weight while maintaining, or even enhancing, structural 
preformance through the use of fiber-reinforced composite materials have been underway for 
approximately 50 years. Although many advancements have been made, and several large-scale 
structures have been built and tested, advances in understanding the mechanics of composite 
structures, especially those with extreme orthotropies and anisotropies, need to continue. One 
common, fundamental structural element is the circular cylinder. For example, understanding the 
buckling response of stiffened cylinders is fundamental to the design of aircraft fuselages and 
launch-vehicle fuel tanks and interstages. Likewise, understanding how edge effects associated 
with end fittings on the tubular members of a lunar lander affect the structural performance is 
significant. When dealing with generally laminated cylinders, deformation modes associated with 
bending, twisting, and stretching can be highly coupled, and failure to fully understand these 
effects on a design can have dangerous and costly unintended consequences. Thus, continued 
advancement of understanding the basic behavioral characteristics of generally laminated circular 
cylinders is an important part of the quest to develop high performance, or special purpose, 
structures by exploiting orthotropies and anisotropies. Moreover, advancing this understanding 
can be beneficial to postmortem analysis of failures that always occur when the limits of design 
technology are pushed, and should help to mitigate risks. Furthermore, advancing this 
understanding may be beneficial to the development of material-characterization test 
methodology for future exotic, advanced materials. Improved understanding of fundamental 
structural elements also provides insight into finite element modeling of complex flight-vehicle 
structures. 

Several analyses for predicting the stresses, strains, and displacements of laminated- 
composite cylinders that are relevant to the present study can be found in the literature. Perhaps, 
one of the first significant papers to address anisotropic cylinders was given by Dong et al. 1 In 
this paper, equations for the linear static response of general shells are given that are based on the 
simplifications attributed to Donnell, Mushtari, and Vlasov (see Brush and Almroth 2 ). In addition, 
the general-shell equations are specialized to circular cylinders and a stress-function formulation 
is given. Moreover, a solution is presented for the axisymmetric bending of a semi-infinite 
cylinder subjected to edge moments. Similarly, an analysis for anisotropic cylinders undergoing 
axisymmetric deformations was given by Gulati and Essenburg 3 in 1967 that is based on a shell 
theory more refined than Donnell-type theories and that includes the effects of transverse shear 
deformation. 

In 1970, Stavsky and Smolash 4 presented a detailed thermoelastic analysis of orthotropic 
composite cylinders. In particular, numerous results are presented for fixed-end cylinders 
undergoing axisymmetric deformations. Also in that year, Vicario and Rizzo 5 presented an exact 
solution for the axisymmetric response of symmetrically laminated cylinders with coupling 
between extension and shear (unbalanced). Later, in 1972, Reuter 6 presented an analytical 
solution for axisymmetric deformation of unbalanced symmetrically laminated circular 
cylindrical shells and balanced unsymmetrically laminated circular cylindrical shells that are 
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subjected to internal pressure and thermal loads. Reuter’s solution is based on Donnell-type linear 
equations and numerical results are presented for filament- wound cylinders made of heat-treated 
carbon-carbon material. 

In 1975, Hennemann 7 presented an analytical solution for the axisymmetric response of 
homogeneous orthotropic cylinders and unsymmetrically laminated cross-ply cylinders subjected 
to uniform axial-compression and lateral pressure loads. Similar results were presented by Jones 
and Henemann 8 later in 1980. In 1984, Cheng and He 9 derived a pair of complex conjugate, fourth- 
order equations, that are based on Flugge’s 10 corresponding equations, for specially orthotropic 
cylinders that can be solved in closed form. Moreover, eigenfunction solutions were presented 
that include the solution for the axisymmetric bending boundary layer and several simplifed 
equations were presented and their relative accuracy was analyzed. 

Chaudhuri et al." presented an analytical solution for unbalanced unsymmetrically laminated 
circular cylindrical shells that are subjected to internal pressure in 1986. Their analysis is based 
on a variant of the Love-Kirchhoff shell theory, 12 which uses an expression for the change in 
surface twist that was given by Timoshenko and Gere. 13 Numerical results are also presented for 
a two-ply shell that demonstrate the coupling effects of the shell anisotropies. Later, in 1989, Abu- 
Arja and Chaudhuri 14 presented similar closed-form solutions for unbalanced, symmetrically 
laminated cylinders and balanced unsymmetrically laminated angle-ply cylinders subjected to 
pressure loads that include the effects of transverse-shear deformation. In 1989, Butler and Hyer 15 
presented an axisymmetric-shell analysis for a general unbalanced unsymmetrically laminated 
circular cylindrical shell that is subjected to axial compression, torsion, or thermal loading. 
Several results are presented that demonstrate the coupling effects of the shell anisotropies on the 
structural response. 

In the early 1990s, Hyer and Paraska 16 17 presented studies of the axisymmeric response of 
balanced symmetrically laminated and balanced unsymmetrically laminated cylinders subjected 
to axial compression and thermal loads. These studies are in the context of nonlinear prebuckling 
deformations that occur as a result of the compression and thermal loads considered. In particular, 
the effects of laminate stacking sequence on the extent and character of the bending boundary 
layers are presented for two groups of three similar laminates. Two of the laminates are 
unsymmetric. Birman 18 presented equations for axisymmetric thermal-stress analysis of ring- 
stiffened orthotropic cylinders in 1991 that include the first-order effects of transverse-shear 
deformation. In 1992, Birman 19 extended this work to symmetrically laminated cross-ply and 
specially orthotropic cylinders, and developed axisymmetric solutions based on a refined shear- 
deformation shell theory. Additionally, Birman 20 extended his work to include nonlinear 
axisymmetric deformations of generally laminated cylinders in 1993. 

In 1997, Vinson 21 presented a linear analysis that focuses mainly on balanced, symmetrically 
laminated cylinders, and an expression is given for the length of the bending boundary layers near 
the cylinder ends that is based on the Love-Kirchhoff shell equations. In contrast, McDevitt and 
Simmonds 22 used Golden veizer’s static-geometric duality principle 23 to reduce the Sanders-Koiter 
equations 24,25 for fully anisotropic, right-circular cylindrical shells to two coupled fourth-order 
equations that use a stress and a curvature function as the unknown, primary field variables. The 
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reduction is accomplished by adding certain negligibly small terms to the stress-strain relations, 
which are intrinsically in error because they must be established experimentally. The approach 
demonstrates how the static-geometric duality principle can be used to reduce greatly the amount 
of algebra needed to obtain results. Eigenfunction solutions are also presented for specially 
orthotropic cylinders that are in agreement with corresponding results presented by Cheng and 
He. 9 Moreover, asymptotic formulas that can be used to determine axisymmetric bending 
boundary layer attenuation lengths and the decay of other non-axisymmetric, self-equilibrated 
edges loads are given. 

An in-depth analytical, parametric study of the attenuation of bending boundary layers in 
symmetrically and unsymmetrically laminated thin cylindrical shells, including balanced and 
unbalanced constructions, was presented for nine contemporary material systems by Nemeth and 
Smeltzer 6 in 2000. Herein, and in reference 26, the term "balanced" refers to laminates that do not 
exhibit coupling between extension or contraction and shearing. The analysis is based on the 
linear Sanders-Koiter shell equations and specializations to the Love-Kirchhoff shell equations 
and Donnell’s equations are included. Two nondimensional parameters are indentified that 
characterize and quantify the effects of laminate orthotropy and laminate anisotropy on the 
bending boundary-layer decay length in a very general and encompassing manner. In this 
particular study, herein, and in much of the recent literature; the term "orthotropy" refers to 
orthogonally anisotropic laminates that exhibit deformation and stress states in which dilatation 
is completely uncoupled from distortion, with respect to the given coordinate system. In contrast, 
the term "anisotropy" implies the presence of material-induced coupling between dilatation and 
distortion. By using these nondimensional parameters, a substantial number of structural design 
technology results are presented in reference 26 that apply to a wide range of laminated-composite 
cylinders. Also in 2000, Paris and Costello 27 presented equations for analyzing cord-reinforced 
composite cylindrical shells and closed-form solutions are given for the axisymmetric response. 
Their analysis is based on Flugge’s equations. 

In 2001, Ye 28 ' 29 presented a three-dimensional elasticity formulation for the decay rates of 
stresses and displacements in cross-ply and angle-ply laminated composite cylinders subjected to 
self-equilibrated end loads. Analyses for the axisymmetric response of semi-infinite laminated 
cylinders that are based on three-dimensional elasticity theory, on a refined shear-deformation 
shell theory, and on classical Love-Kirchhoff shell theory were presented by Bhaskar and 
Balasubramanyam 30 in 2002. Somewhat unique to their study is a sizeable collection of tabular 
results in addition to figures. In 2003, Paris and Costello 31 extended their previous work on cord- 
reinforced composite cylindrical shells to include generally laminated wall constructions. This 
later analysis is also based on Flugge’s equations and closed-form solutions are given for the 
axisymmetric response. Also in 2003, Kollar and Springer 32 gave a relatively terse presentation 
of the general closed-form solution for the axisymmetric response of fully anisotropic cylinders 
subjected to mechanical and thermal loads that are also based on Flugge’s equations. 

The literature described previously herein contains many elements of a general formulation 
for the axisymmetric behavior of cylindrical shells. However, none of the studies is complete and 
many details of the solutions are omitted. The objective of the present study is to give a detailed, 
complete formulation for axisymmetric deformations of a generally laminated circular cylindrical 
shell subjected to combined mechanical and thermal loads. As a first element of a building block 
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process to obtain improved theories, the analysis presented herein is based on the linear Sanders- 
Koiter shell equations, which neglects transverse-shear deformations. Moreover, the analysis is 
an extension of that presented by Nemeth and Smeltzer 26 in that it uses positive-definiteness of the 
strain energy density to determine the specific form of the analytical solution. This approach 
results in a significant simplification of the analysis procedure. To accomplish this objective, the 
Sanders-Koiter equations for axisymmetric deformations are given and then manipulated into a 
convenient form that facilitates obtaining the analytical solution. Then, a systematic presentation 
of the exact solution is given, followed by examples that demonstrate the typical solution 
processes. Approximate solutions are also presented to add credibility to the present study and to 
illustrate potential pitfalls of using approximate solutions in the analysis and design of laminated- 
composite cylinders. 


Equations Governing Axisymmetric Behavior 

The ordinary differential equation that governs the axisymmetric bending behavior of a right- 
circular cylinder that is subjected to edge loads or displacements and surface tractions is obtained 
by first specializing the linear Sanders-Koiter shell equations, that are given in Appendix A, for 
axial symmetry. For the equations presented herein, x and 0 denote the axial and 
circumferential coordinates of a right-circular cylinder, respectively, as shown in figure 3. The z- 
coordinate shown in the figure is positive in the radially outward direction, and z = 0 corresponds 
to the middle surface of the shell. The cylinder has length L, radius R, and uniform thickness h. 
The specialization to axial symmetry is obtained by eliminating all terms in the equations that are 
differentiated with respect to the circumferential coordinate, 0. The resulting set of equations for 
axisymmetric behavior is given in Appendix B. 

The particular equations used in the present study to obtain the exact solution for 
axisymmetric behavior are obtained by first noting that integration of equation (Bl) yields 


N x (x) = N 0 + N,(x) (la) 

where N 0 is a constant of integration and 

N,(x) = - J q x (x) dx (lb) 

Next, equations (B2) and (B5) are combined to get 

N x0 '(x) + -g/cj + ^c 2 |m xB '(x) + q 0 (x) = 0 (2) 

where the prime mark denotes differentiation with respect to the independent variable, x. For 
convenience, the parameter 
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( 3 ) 


M- — c ] + 2^2 

is introduced such that the Sander s-Koiter equations are given by q = 3/2 and the Love-Kirchhoff 
equations are given by u = 1 . Donnell’s equations are given by q = 0. Similarly, the function 

T(x) = N x0 (x) + ^M x0 (x) (4) 

is introduced so that equation (2) becomes 


T'(x) + q e (x) =0 


( 5 ) 


and the corresponding boundary condition given by equation (B19) becomes 


T = T or v = v 


( 6 ) 


at x = 0 and x = L. Integration of equation (5) yields 


f(x) = T 0 + T,(x) 


(7a) 


where T 0 is a constant of integration and 


F ' (X) = "J 


q e (x) dx 


(7b) 


Next, equations (B3) and (B4) are combined to get 


M, (xl - 


Nix 


R 


+ q„(x) = o 


( 8 ) 


The next step in the analysis is the simplification of the constitutive equations. First, by using 
equations (B8) and (3), equation (B14) is expressed as 


o _ q o 

K x0 - *e 


( 9 ) 


By using equations (4) and (9), the constitutive equations given by matrix equation (B15) are 
expressed as 
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(12b) 

h + ^Ir) h 2 
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h ) D 16 
Ri h 
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.R M Th 
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where the subscripted A, B, and D terms appearing in these equations are the usual laminate 
stiffnesses of classical plate and shell theory (see reference 34), in which the fiber orientation of 
a typical ply is denoted by the angle ()) shown in figure 3. The motivation for writing the 
constitutive equations in this form is that the matrix equation given by equation (10) is the only 
part of the full constitutive equations that appears in the strain energy density function, which is 
used in the present paper to determine the corresponding positive-definiteness conditions. Noting 

that e° = u r (x) and 'ha - v r (x) are the only strains in equation (10) that depend on the axial 
displacement u(x) and the circumferential displacement v(x), it is useful to extract 
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(13a) 


from equation (10). Using equations (la), (7a), (B6) - (B8), and (B12) with equation (13a) gives 


/n o + N,(x) + N™(x)| _ A n A 16 

\t 0 + T,(x) + T Th (x) / [ S ,6A 66 
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(13b) 


Solving for the derivatives of the axial displacement u(x) and the circumferential displacement 
v(x) gives 
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(15b) 

(15c) 


and where w(x) is the radial deflection that is positive-valued when outward. Equations (14) and 
(15) indicate that the circumferential displacement v(x) becomes uncoupled from the axial 

displacement u(x) and the radial displacement w(x) when A 16 = A 26 = B 16 = 0 , which implies 
that A l6 = A 26 = B 16 = D i6 = 0. In addition, the constitutive equation, equation (10), indicates that 

N x , N 0 , and M x become uncoupled from the torsional, shear strain Y°e when A 16 = A 26 = B 16 = 0 , 

and that T , that is defined by equation (4), becomes uncoupled from e°, and k°. Furthermore, 

equation (11) indicates that M e becomes uncoupled from Y°e when B 26 = 0, which implies B 26 = 
D 26 = 0. 


Next, equations (14) are then substituted into e ° = u, ( x ) and Yxe = v r ( x ) , and the resulting 

expressions for £" and Y°e , along with equations (B7) and (B12) are substituted into the 
constitutive equations given by equation (10). This action converts the strains and stress 
resultants in equation (10) into functions of the radial displacement w(x). In particular, 

N e (x) = A 12 u'(x) + A 22 ^> + A 26 v'(x) - B 12 w"(x) - n”(x) (16a) 


becomes 


N e (x) - (A 12 5 u + A 26 5 16 )[N 0 + N,(x) + N x (x)] + (A 12 * 16 + A 26 5 66 )[t o + T,(x) + T (x)] 


+ ( A 22 + A i2 d i2 + A 26 ^ 62 j p (B 12 + A 12 e 12 + A 26 6 62 Jw (x) N e (x) 


(16b) 


R 


and 


M x (x) = B,,u'(x) + + B IC v (x) 


W M , 5 -O 

’ R 


D„w"(x) 


M v 


(17a) 
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becomes 


M x (x) - (B n « u + B 16 5 16 )[N 0 + N,(x) + N x (x)] + (B u * 16 + B 16 5 66 )[t o + T,(x) + T (x)] 
+ (B IZ + B n * I2 + B 16 ^ 62 ) (D u+ B n * 12 + B 16 * 62 )w (x) — M x (3 


Substituting equations (16b) and (17b) into equation (8) yields 

C, w""(x) + C 2 w"(x) + C 3 w(x) = C 4 (x) 
The constant coefficients are given by 


( 18 ) 


R' 


and e = 1 + B 1^12 + B 16«62 
D„ 

(19) 

12 A 12^ 12 26^ 62 ^11^12 ^16^62) 

(20) 

2 ( A 22 + A l2 ^ 12 + AJ 62 j 

(21) 


R 


The function C 4 (x) is given by 


C 4 (x) = q„(x) - M x '’"(x) + ^ Ng h (x) 

- i(A 12 * n + A 26 5 16 )[N 0 + Nj(x) + N x h (x) j - A 

+ (b„5h + Bjgijgj^N/^x) + N x "(x)] + (B m 5 16 + B 16 5 66 j 


12*16 " I " A 26 5 66 


T 0 + T,(x) + T (x) 
T/'(x)+f Th "(x)] 


( 22 ) 


Equations (19) - (21) are expressed in terms of elements appearing in matrix equation (10) further 
by using equations (15); that is, 


C = - — 
2 R 


B 12 + 


(A| 6 A 26 A 12 A 66 )B u + (A 12 A 1s AjjA^jB] 

a„a 66 -a 16 


(23) 


C = — 

^3 2 

R 


AhA 22 A 12 ]a 66 A u A 26 A 22 A 16 + 2 A 12 A 16 A 26 


A n A 66 -A 16 


(24) 
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(25) 


t A 66 Bi, + A 11 b; 6 -2A„B, i B i 


A hA 66 A 16 Dj 


Next, equation (18) is divided by the coefficient C, given by equation (19) to get 

w""(x) + 4S w"(x) + 4Q w(x) = P(x) 
where the constants S and Q are given by 


S = 


1 


4C, 


2RD,,e 


B, 2 + 


(A 16 A 26 A 12 A 66 )B u + (A 12 A 16 A u A 26 )Bj 

a„a 66 -a 2 16 


Q = 


"3 V 


A u A 22 A 12 A 66 A u A 26 A 22 A 16 + 2 A 12 A 16 A 26 


4C, 


4R D u c I A u A 66 — A lf 


The function P(x) is given by 


C 4 (x) 

P(x) = ^ Po + P| (x) 

I 


where 


Pn = - 


[(A, 2 «h + A 26 a 16 )N 0 + (A 12 5 16 + A 26 « 66 )T 0 ] 


rd,i e 


Pi(x) = — - — [q„(x) - M™"(x) + ± N^x)] 

Du e 

[( A 12 « n + A 26 « 16 )(n i (x) + N x (x)) + (A 12 * 16 + A 26 5 66 )(T[(x) + T (x))j 

RD„ e 

+ " (B„*n + B 16 * 16 )[n"(x) + N^'W] + -J— (B n i I6 + B 16 * 66 )[t/'(x) + f Fh "( 

Du & D u ^ 


(26) 

(27) 

(28) 

(29a) 

(29b) 

(29c) 
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Exact Solution 


Equation (26) is a linear, nonhomogeneous, fourth-order ordinary differential equation with 
constant coefficients, and the form of its solution depends on the numerical values of the constants 
S and Q. To determine the specific form of the solution, it is useful to examine the positive- 
definiteness conditions on the strain-energy density function. The strain energy density function 
71 for this problem is given by (e.g., see Paraska, 17 p. 16) 


2 71 = (N, - N™)e° + (N B - NiTK + (N xe - + (M x 


, ,Th\ o 

M x )k x 


+ 


(M x 


a * Th \ o 

M x 0 )k x 


(30) 


By using equations (4) and (9), the strain-energy density function is expressed as 


2 71 = (N x - NTK + (N e - + (f - T T, ’) Y : e + (M x - mT)< 


(31) 


The strain energy density is expressed in terms of the total strains, mechanical plus thermal, and 
constitutive terms by using the constitutive equation given by equation (10); that is. 


7t = 


1 

2 



A ii A 12 A 16 B n 

A A A R 

n 12 n 22 -^26 u 12 


A 16 A 26 A 66 B 16 


B,, B 12 B 16 D„ 



(32) 


The stiffness terms in equation (32) that have overbars are defined by equations (12) and are 
functions of the shell wall thickness-to-radius parameter, h/R. Because the strain energy density 
function is positive valued, the matrix appearing in equation (32) is positive definite. As a result 

(e.g., see Zwillinger 33 , pp. 133-134), the diagonal terms A n , A 22 , A 66 , andD n are positive-valued, 
and the following determinants are positive valued 

= A n A 22 - A 2 12 >0 (33) 


A ii A 12 

A A 

2 *- 12 2^-22 


A ii A 12 A 16 

A 12 A 22 A 26 
A 16 A 26 A 66 


— (a u A 2 2 a l2 ) A 66 A|iA 2 6 A 22 A 16 + 2 A i 2 A 16 A 26 > 0 (34) 


in addition to the determinant of the constitutive matrix in equation (32). Moreover, by 
rearranging the strain energy density function into the form 
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( 35 ) 





11 A i6 A 12 B n 
A i6 A 66 A 26 B 16 

A A A R 

• rv 12 ■ rl 26 ^ 22 JJ 12 

B„ B 16 B 12 D n 



the following additional positive-definiteness condition is obtained 


A„ A 16 
a i 6 a 66 


= a 11 a 66 - a 16 >o 


(36) 


The homogeneous solution for equation (26) involves the square root of the quantity Q - S 2 . 
By using equations (27) and (28), this quantity is given by 

Q - S 2 = 4ClC \ C2 (37) 

16C, 

Substituting equations (19) and (23) - (25) into equation (37) and simplifying, the quantity Q - S 2 
is found to be given by 


Q-S = 


4C'r 2 (a u A 66 - A 2 16 ) 


An A 12 A 16 B n 

A 12 A 22 A 26 B 12 
A 16 A 26 A 66 Bjg 

Bn B j 2 B ,6 D n 


(38) 


It follows logically, that Q - S 2 > 0 because the positive definiteness of the strain energy density 
function requires that the determinants in equations (36) and (38) be positive valued. Moreover, 
Q - S 2 > 0 implies that Q > 0, and Q > 0 implies that C 3 /C, > 0. Equations (34), (36), and (24) 
indicate that C 3 >0. Thus,C 3 /Cj>0 yields the condition that C, = D n e>0 (see equations (19) 


and (25)). Because D u >0,e> 0. To reinforce the postive valuedness of Q, it is convenient to 
introduce the expression 


T = Q = 


AnA 2 2 A 12 I A 66 A,jA 26 A 22 A 16 + 2A i 2 A i 6 A 2 



( 2 \ 

4R d,, e 

a„a 66 - a, 6 


l ) 


(39) 


such that T 2 - S 2 > 0, and to express equation (26) as 

w ;,,, (x) + 4S w"(x) + 4T~ w(x) = P 0 + P,(x) (40) 
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Equation (40) is a linear, fourth-order, nonhomogeneous ordinary differential equation with 
constant coefficients. The characteristic equation of equation (40) is given by 

X 4 + 4S X 2 + 4T 2 = 0 (41) 

Using the knowledge that T 2 - S 2 > 0, the roots of the characteristic equation are obtained from the 
quadratic formula; that is, 

(X 2 ),, = 2(— S ±i V T 2 -S 2 ') (42) 

where i = . Solution of this equation for X yields four roots of equation (41) that are pairs 

of complex conjugates that are given by 


x 1AM = ± ( U T - 8 * * V T + S ) 


(43) 


The homogeneous solution of equation (40) can be written as follows 


where 


K/,( 

x) + K/,(x) + K 3 ^ 3 (x) + K/ 4 (x) 

(44) 

^,( x 

) = e sin( VT + S x) 

(45a) 


) = e cos( 4T + S x) 

(45b) 


I = e * sin( VT + S x) 

(45c) 

^ 4 ( X S 

1 = e ’ cos( ZrT S x) 

(45 d) 


are basis function spanning the solution space and where 0 < x < L. The symbols K p K,, K 3 , and 
K 4 are real-valued constants that are determined from the boundary conditions. The solution 
given by equation (44) represents a damped, oscillatory response that decays from each end of the 
cylinder. The regions near the edges of the cylinder, where the amplitude of w H (x) typically varies 
the most, are called the bending boundary layers. This particular form for the basis functions was 
selected to prevent numerical issues associated with evaluating exponential functions with a large 
positive-valued argument. 

Next, it is important to observe that the constant-valued term P 0 in equation (29a) contains 
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two of the unknown constants of integration, N 0 and T 0 . Thus, it is convenient to partition the 
particular solution of equation (40) into 


w p (x) = w" + w p (x) 


(46) 


where w° is a constant obtained as the particular solution associated with P 0 and w p (x) is the part 
of the particular solution associated with P , ( x ) . Moreover, it is convenient to express 


^ = R(crN 0 + c t T 0 ) 


(47a) 


where 


A «A 26 A 12 A 66 


A|iA 2 2 A 12 ]a 66 A u A 26 A 22 A 16 + 2A 12 A 16 A, 


(47b) 


c T - 


A|eA 12 A u A 2 


A n A 22 A | 


A 


(47c) 


A 66 A u A 26 A 22 A 16 + 2A 12 A 16 A 2 


The form of w p (x) depends on the specific functional representations of the applied surface 
tractions and thermal fields. By using equations (44), (46), and (47a), the general solution to 
equation (40) is expressed as 

w(x) = K/,(x) + K/ 2 (x) + K 3 ^ 3 (x) + K/ 4 (x) + R(c r N 0 + c t T 0 ) + w p (x) (48) 

With the form of w(x) known, the first of equations (14) is expressed as 

u'(x) = U 0 '(x) + U/(x) K, + U 2 '(x) K 2 + U/(x) K, + U 4 '(x) K 4 + U 5 'N 0 + U/T 0 (49) 


where 


U 0 '(x) = 5 h [n ,( x) + N™(x)] + 5 16 [t,(x) + T Th (x)] + - 5 12 w p "(x) (50a) 

- ^k( x ) 

U k ( x ) = ^ 12 — ^ «i 2 ^ k ( x ) with k = 1, 2, 3, and 4 (50b) 

U 5 = «u + 4 i 2 Cfj (50c) 
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U 6 — *16 + ^ 12^T 


(50d) 


Likewise, the second of equations (14) is expressed as 

v'(x) = V 0 '(x) + V/(x) K, + V 2 '(x) K 2 + v 3 '(x) K 3 + V 4 '(x) K 4 + V 5 'N 0 + V 6 ' 

where 

V 0 '(x) = « 16 [n,(x) + NlV)] + * 66 [t i (x) + T (x)] +^ 62 -^-^ -« 62 W p "(x) 

^ (x) 

V k '(x) =^ 62 — ^ tf 6 2 ^/ , (x) with k = 1, 2, 3, and 4 

V/ = * 16 + ^ 62 c N 

V 6 — *66 ^62^T 

Next, the derivatives of the basis functions, given by equations (45), are expressed as 

*/(x)= vTTs/ 2 (x)- vT^sY,(x) 

* 2 '(x) = - vTTS^x)- vT^S ^ 2 (x) 

* 3 '(x)= vTTs^ 4 (x) + vT^S^ 3 (x) 

#/(x) = - vTTsY 3 (x) + VT^S^ 4 (x) 

4 /'(x) = - 2S ^ ,(x) - 2 \/T-S 2 ' 4 2 (x) 

<T 2 "(x) = 2 Vt-S 2 ' *,(x) - 2S ^ 2 (x) 

<£ 3 "(x) = - 2S * 3 (x) + 2 Vt 2 - S 2 * 4 (x) 

^ 4 "(x) = - 2 H - S 2 ' * 3 (x) - 2S * 4 (x) 


o (51) 

(52a) 

(52b) 

(52c) 

(52d) 

(53a) 

(53b) 

(53c) 

(53d) 

(53e) 

(53f) 

(53g) 

(53h) 
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By using these expressions for the derivatives, equations (50b) are expressed as 



^ + 2S5 I2 


,ix) + 2a 12 Vr-S 2 ' * 2 (x) 


(54a) 


U 2 '(x) = -25 12 Vt 2 -S 2 ' *,(x) + 



+ 2Sa 12 


^ 2 (x) 


(54b) 


U/(x) = 


± + 2Si n 


^ 3 (x)-25 i2 Vt 2 -S 2 '^ 4 (> 


(54c) 


U/(x) = 2c 12 Vt-S 2 ' <£ 3 (x) + 



+ 2S c, 


(54d) 


Likewise, equations (52b) are expressed as 


V/(x) = 


^ + 2S5 62 


’ 1 (x) + 25 62 v / T 2 -S 2 ' ^ 


(55a) 


V/(x) = -25 62 Vt 2 -S 2 ' *,(x) + 



+ 2S« 62 


*,(x) 


(55b) 


V 3 '(x) - 


-jf + 2 Sc 6 


^ 3 (x) - 25 62 Vt 2 -S 2 ' 


(55c) 


V/(x) = 25 62 n/t 2 -S 2 ' * 3 (x) + 



+ 2S« fi 


(55d) 


Next, equations (49) and (51) are integrated directly to get 


u( x ) = U 0 (x) + U,(x) K, + U 2 (: 


x) K 2 + U 3 (x) 


w 
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and 


where 


u 0 (x) = 5„ J [n,(x) + Nr ( x)]dx + s I6 J [Tj(x) + T Th (x)]dx + ^ J 


l=r I w (x) dx - c 12 w '(x) 


V 0 (x) 


: «i6 J [n,(x) + N™(x)]dx + <* 66 J" [t,(x) + T (x)]dx + ^J* w p 


-r 1 I w p (x) dx - i 62 w p (x) 


U k (x) 

V k (x) 



U k '(x) dx 
V/(x) dx 


(58a) 


(58b) 


(58c) 


(58d) 


for values of k = 1, 2, 3, and 4; and u and v are the corresponding constants of integration. 
Next, the integrals of the basis functions, given by equations (45), are expressed as 

J ^(x)dx = - ^ r (yT^SY,(x)+ vTTS^x)) (59a) 

| * 2 (x)dx = ^(-^r = SV 2 (x) + VTTS *,(x)) (59b) 

J * 3 (x) dx = Jf( * 3 (x) - VTTS * 4 (x)) (59c) 

^ d 4 [x) dx = Jp( ^ 4 (x) + VTTS * 3 (x)) (59d) 

By using these expressions for the integrals, equations (58c) and (58d) are expressed as 

U 1 (x) = c^ 1 (x)-c“/ 2 (x) (60a) 

U 2 (x) = c"^ 1 (x) + c"^ 2 (x) (60b) 
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U 3 (x) = - c u / 3 (x) - cj 4 (x) 


(60c) 


with 


and 


with 


U 4 (x) = c/ 3 (x)-c 1 , 


(60d) 


c“ = 


2T 


25 12 T - 


(60e) 


^2 _ 


/TTS 

2T 


2«i 2 T + ^ 


(60f) 


V 1 (x) = c^,(x)-c^ 2 (x) 

(61a) 

v 2 (x) = c^ 1 (x) + c^ 2 (x) 

(61b) 

V 3 (x) = -c^,(x)-c 2 ^ 4 (x) 

(61c) 

V 4 (x) = c^ 3 (x) -cX 4 (x) 

(6 Id) 


c, = 


2T 



(61e) 


^2 _ 


vTT" S 
2T 


2tf 62 T + 



(6 If) 


Equations (48), (56), and (57) define completely the displacement fields in terms of eight 
constants that are determined from the boundary conditions. These constants are 

{Kj, K,, K 3 , K 4 , u, v, N 0 , T 0 } (62) 


To determine the unknown constants given by equation (62), the boundary conditions stated 
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by equations (B18) - (B21) must be enforced. To enforce these equations, expressions for the 

° 

rotation Ih and the stress resultants are needed. First, the expression for the slope is obtained by 
differentiating equation (48) and using the expressions for the derivatives of the basis functions 


given by equations (53). This task yields 

P°( x ) = |3 0 ( x ) + K,p,(x) + K 2 (3 2 (x) + K 3 p 3 (x) + K 4 p 4 (x) (63) 

where 

Po(x) = - w/(x) (64a) 

p 1 (x)=/r^l^ l (x)-/r+¥^ 2 (x) (64b) 

p 2 (x)=/mT(x)+/r^T(x) (64c) 

|3 3 (x) = - VT^T? ^ 4 (x) - vT^ ^ 3 (x) (64d) 

p 4 (x) = vT+ir <? 3 (x) - <? 4 (x) (64e) 


Next, equations (48), (49), and (51) are substituted into equation (17a) to get 

M x = M x0 (x) + M xl (x) K, + M x2 (x) K 2 + M x3 (x) K 3 + M x4 (x) K 4 + M x5 N„ + M x6 T 0 (65) 


where 


M x0 (x) = B n U 0 '(x) + B 


wJx) _ , 

P +B, f V/( 


R 


d u w; 


M x (x) 


(66a) 


g 

M Ik (x) = BjT/fx) + Y^ k (x) + B 16 V k '(x) - Dn^/'jx) w ph k=l,2, 3, and 4 (66b) 
M x5 = B n U 5 ' +B 12 c R + B 16 V 5 ' (66 c) 

M x6 = B „U S ' + B 12 c T + B 16 V 6 ' (66d) 

Using equation (50a) for U 0 f (x), equation (52a) for V 0 , (x), equations (54) for U/(x) through U 4 ' '(X), 
equations (55) for V/(x) through V/(x), equation (50c) for U/, equation (50d) for U 6 \ equation 
(52c) for V/, and equation (52d) for V/ , and equations (53) for the second derivatives of the basis 
functions gives 
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(67a) 


M x0 (x) — (BjjSji + B 16 5| 6 )[n,(x) + N x (x)j + ^B n 5| 6 + B 
+ (B 12 + Bj u + B,/ 62 )^^ (D u + B n tf 


i6*66)[t,(x) + T (x)] 

12 + B l 6 tf 62 )w p "(x) - M' h (x) 


where 


M xl (x) = c^V,(x) + c^V 2 (x) 

(67b) 

M s 2 (x) = -c^ 1 (x) + c^ 2 (x) 

(67c) 

M x 3 (x) = c^ 3 (x) -cf^Jx) 

(67d) 

M x4 (x) = c“^ 3 (x) + c^ 4 (x) 

(67e) 

!M x5 = c r (B 12 + B n ^ 12 + B,/ 62 ) + Bjjctjj + B 16 * 16 

(67f) 

-M x6 = Ct(B 12 + Bj l2 + B 16 ? 62 ) + B„, 16 + B 16 <H 66 

(67g) 

M B I / I2 + B 12 + B 16 ^ 62 / _ — x 

Cj — ^ + 2S^B n c 12 + B 16 £ 62 + Dn J 

c“ = 2 Vt 2 - S 2 (B„5 I2 + B I6 e 62 + D u ) 

(67h) 

(67i) 


Next, the transverse shear stress resultant Q x (x) is obtained by using equation (65) with equation 
(B4) to get 

Q x (x) = M/(x) . Q x0 (x) + Q xl (x) K, + Q x 2 (x) K 2 + Q x 3 (x) K 3 + Q x4 (x) K 4 (68) 

where Qxk( x ) = M xl /(x) for values of k = 1, 2, 3, and 4. Differentiating equation (67a) gives 

QxoW = M x0 '(x) = (B ll5ll + B 16 5 16 )[n/(x) + N X V)] + (B u 5 16 + B 16 ^ 66 )[t/(x) + T "'(x)] 

_ ,, , (69a) 

+ [B 12 + Bj I2 + BJ 62 ] — ^ [D u + B„tf 12 + B^Jw/^fx) - M x '(x) 

Similarly, differentiating equations (67b) - (67e) and using equations (53) gives 
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Q xl (x) = - [c“ VT^S + c 2 “ vTTs Y , ( x ) + [c“ /FT S - c“ vT^S ]^ 2 (x) (69b) 

Q x2 (x) = [c“ vT^S - c“ vTTs ,(x) - [T/FTs +cr/r^s]^ 2 (x) (69c) 

Q x3 (x) = [ c “ vT^S + c“ /FTs]^ 3 (x) + [cr /FT S' - c“ vT^s]^ 4 (x) (69d) 

Q x4 (x) = [cr /T-S - cr /FTs ]^ 3 (x) + [cr /FTs + cr /F^S ]^ 4 (x) (69e) 


All expressions needed to enforce 8 of the 16 boundary conditions given by equations (B18)- 
(B21) have been obtained. In particular, equations (la) and (56) are used to enforce the boundary 
conditions stated by equation (B18). Equations (7a) and (57) are used to enforce the boundary 
conditions stated by equation (6), and equivalently equation (B19). Equations (68) and (48) are 
used to enforce the boundary conditions stated by equation (B20), and equations (63) and (65) are 
used to enforce the boundary conditions stated by equation (B21). The resulting 16 possible 
boundary conditions, in which one must be selected from each pair, are listed as follows: 


u( 0 ) = u 0 : u,( 0 ) K, + u 2 (0) K 2 + u 3 (0) K 3 + u 4 (0) K 4 + u = u 0 - U 0 (0) or (70a) 

N,(0) = N 0 : N 0 = N 0 - N,(0) (70b) 

v(0) = v 0 : V i (0)K 1 + V 2 (0)K 2 + V 3 (0)K 3 + V 4 (0)K 4 + v = v o -V o ( 0) or (71a) 

T(0) = f„: T 0 = f 0 - Tj(o) (71b) 

w(o) = w 0 : k/ 1 { 0 ) + k^ 2 ( 0 ) + k/ 3 ( 0 ) + k/ 4 ( 0 ) + r(c s n 0 + c t T 0 ) = w„-w p ( 0 ) 0 r (72a) 
Qx(0) = V o : Q xl (0) K, + Q x2 (0) k 2 + Q x 3 (0) K, + Q x4 (0) K 4 = V 0 - Q x0 (0) (72b) 

P°(o) = p 0 : K 1 p I (0) + K 2 p 2 (0) + K 3 p3(0) + K 4 p 4 (0) = p 0 -| 3 0 ( 0 ) or (73a) 
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(73b) 


M x (0) = M 0 : 

M xl (0) K, + M x2 (0) K 2 + M x3 (0) K 3 + M x 4 (0) K 4 + M x5 N 0 + M x6 T 0 = M„ - M x0 (0) 


u(L) = u L : 

U,(L) K, + U 2 (L) K 2 + U 3 (L) K 3 + U 4 (L) K 4 + (U 5 'N 0 + U 6 'T 0 )L + u = u L - U 0 (L) 

N x (L) = N l : N 0 = N l - N,(L) 

v(L) = v L : 

V , (L) K, + V 2 (L) K 2 + V 3 (L) K 3 + V 4 (L) K 4 + (V 5 'N 0 + V 6 'T 0 )L + v = v L - V 0 (L) 

T (x) = T L : T 0 = T L - T,(L) 


or (74a) 
(74b) 
or (75a) 
(75b) 


W (L) = w L : 


K/ t (L) + K/ 2 (L) + K/ 3 (L) + K/ 4 (L) + R(c n N 0 + c t T 0 ) = w L - w p (L) 


or (76a) 


Q x ( x ) — V L : Q XI (L) K, + Q x2 ( l ) K 2 + Q x3 (L) K 3 + Q x4 (L) K 4 — V L - Q x0 (L) 


(76b) 


(3° x (L) = p L : K l (5 1 (L) + K 2 (5 2 {L) + K 3 (3 3 (L) + K 4 p 4 (L) = p L -p 0 (L) or (77a) 


M X (L) = M l : 

M x1 (L) K, + M x2 (L) K 2 + M x3 (L) K 3 + M x4 (L) K 4 + M x5 N 0 + M x6 T 0 = M l - M x0 (L) 


(77b) 


For certain combinations of boundary conditions, rigid-body displacements may be present that 
must be eliminated in order to have a unique solution. 

After determining the unknown constants given by equation (62), the displacements, strains, 
and stress resultants are known. The stress resultant N fl is obtained from equation (16a). 
Substituting equations (48), (49), and (51) into equation (16a) gives 


No = Ng 0 (x) + N fll (x) K, + N 02 (x) K 2 + Ng 3 (x) K 3 + N e4 (x) K 4 + N 95 N 0 + N fl6 T 0 (78) 
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where 


(79a) 


Neo(x) = A I 2 U 0 '(x) + A 22 ^— + A 26 V 0 '(x) - B 12 W p "(x) - Ne“(x) 


N flk (x) = A I 2 U k '(x) + ^^ k (x) + A 26 V k '(x)-B 1 / k "(x) with k= 1,2, 3, and 4 (79b) 

n 65 = a, 2 u 5 ' + a 22 c r + A 26 V 5 ' ( 79 c) 

N e6 = A 12 U 6 ' + A 22 c t + a 26 v 6 ' ( 79 d) 

Using equation (50a) for U/(x), equation (52a) for V 0 '(x), equations (54) for U/(x) through U/(x), 
equations (55) for V/(x) through V 4 ' '(X), equation (50c) for U 5 \ equation (50d) for U 6 \ equation 
(52c) for V 5 \ and equation (5 2d) for V/ , and equations (53) for the second derivatives of the basis 
functions gives 


Neo( x ) — (A 12 * u + A 26 *, 6 )[n,(x) + N x (x)] + (A l2 # 16 + A 26 5 66 )[t,(x) + T (x)j 


, w„ x 


(80a) 


+ (A22 + A , 2 ^ | 2 + A 26 ^ 62 ) (B 12 + A l 2 i l2 + A 26 5 62 jw p (x) N e (x) 


N e i(x) = ,(x) + c 4 ^ 2 (x) 

(80b) 

N 02 (x) = - c 4 ^,(x) + c 2 ^ 2 (x) 

(80c) 

N 63 (x) = c^ 3 (x) -c 4 ^ 4 (x) 

(80d) 

N 64 (x) = c 4 V 3 (x) + c^ 4 (x) 

(80e) 

n(A 22 + A 12 £ 12 + A 26 # 62 ) + A 12 * u + A 26 (* 16 

(80f) 

r(A 22 + A i2 i 12 + A 2 / 62 ) + A I2 d l6 + A 26 a 66 

(B0g) 


where 


n A l2 d 12 + A 22 + A 2 / 62 


c, = 


R 


+ 2S(A 12 tfi2 + A 26^ 62 ^ 12 ) 

(81a) 

,C l2 + A 2 6 ^ 6 2 + B 12 ) 

(Bib) 
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Similarly, The stress resultant M e is obtained from equation (1 1) by using equations (B6) - (B8) 
and (B12); that is 


M„ = B I2 u'(x) + B 22 "^ + B 26 v'(x) - D 12 w"(x] - Me h (x) (82) 

Substituting equations (48), (49), and (51) into equation (82) gives 

M a = M 0 O (x) + M ei (x) K, + M 02 (x) K 2 + M 03 (x) K 3 + M e 4 (x) K 4 + M e 5 N 0 + M 06 T O (83) 

where 

M 0O (x) = B 12 U 0 '(x) + B 22 ^l + B 26 V 0 '(x) - D 12 w p "(x) - M 0 h (x) (84a) 

g 

M ek (x) = B 12 U k '(x) + ^^ k (x) + B 26 V k '(x)-D 1 / k "(x) with k=l,2, 3, and 4 (84b) 

Mes = B i2 U 5 ' + B 22 c n + B 26 V 5 ' (84c) 

M e6 = B 12 U/ + B 22 c t + B^v/ (84d) 

Again, using equation (50a) for U/(x), equation (52a) for V 0 , (x), equations (54) for U/(x) through 
U/(x), equations (55) for V/(x) through V 4 '(X), equation (50c) for U 5 \ equation (50d) for U 6 \ 
equation (52c) for V 5 \ and equation (52d) for V/ , and equations (53) for the second derivatives 
of the basis functions gives 

Meo( x ) = (B 12 5 u + B 26 5| 6 )[n,(x) + N x (x)] + (B 12 5 16 + B 26 5 66 )[t!(x) + T (x)] 

w (x) _ rh ( 85a ) 

+ (B 22 + B i2 ^ 12 + B 2( / 62 ) p- (D 12 + B 12 5 12 + B 26 5 62 )w p (x) — M 0 (x) 

M 61 (x) = c“^,(x) + c“^ 2 (x) (85b) 

M 02 (x) = - cy ,(x) + c^V,(x) (85c) 

M 93 (x) = c 3 V ,(x) - c“^ 4 (x) (85d) 

M 04 (x) = cy 3 (x) + C y 4 (x) (85e) 
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= Cf}(] 


^ 26^62 

) + B i 2 5 u + B 26 « i6 

(85f) 

^ 26 ^ 62 ^ 

+ B 12 «| 6 + B 26 « 66 

(B5g) 


where 


M ® 12^ 12 ^22 "I" ® 26"' 62 


C, = 


R 


• / — \ 

+ 2S(B 12 a 12 + B 26 a 62 + D 12 j 

(85h) 

12^12 B 26 c 62 + D 12 ) 

(85i) 


The stress resultant M x0 (x) is obtained by using equation (B15) and (9) to get 

-M-xe( x ) = B l6 £ x + B 26 £ e "T B 66 y x e + D, 6 k x — M x0 

where 

B„ = B„ + n(i)% 

Then, using equations (B6) - (B8) and (B12), equation (86) becomes 

Mix) = B 16 u'(x) + B 26 ^l + B 66 v'(x) - D 16 w"(x) - M x 0 (x) 

Substituting equations (48), (49), and (51) into equation (88) gives 

M x0 = M x 0 O (x) + M x 01 (x) K, + M x 02 (x) K 2 + M x 03 (x) K 3 + M x 04 (x) K 4 + M X 05 N O + M X 06 T„ 

where 

M x 0 O (x) = B 16 U 0 '(x) + B 26 — ^ + B 66 V 0 '(x) - D 16 w p "(x) - M x 0 (x) 

B 


( 86 ) 


(87) 


( 88 ) 


(89) 


(90a) 


M * 8 k ( x ) = B 16 U k '( x ) + ^ k ( x ) + B 66 V k '( x ) - D 16 tf k "( x ) with k=l,2, 3, and 4 (90b) 


M x05 — B 16 U 5 ' + B 26 c s + B 66 V 5 ' 


Mxee — B 16 U 6 ' + B 26 c t + B^V/ 


(90c) 

(90d) 
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Again, using equation (50a) for U/(x), equation (52a) for V 0 , (x), equations (54) for U/(x) through 
U/(x), equations (55) for V/(x) through V 4 '(x), equation (50c) for U 5 \ equation (50d) for U 6 \ 
equation (52c) for V 5 \ and equation (52d) for V/ , and equations (53) for the second derivatives 
of the basis functions gives 


M xeo (x) - (B 16 S u + B 66 5, 6 )[Ni(x) + N x (x)] + (B 16 * i6 + B 66 S 66 )[t,(x) + T (x)] 


(B 26 + B 


16^ 12 "t" B 6 / 62 ) 


R 


(9 16 B 16 tf 12 + B 66 £ 62 Jw p (x) M x0 (x) 


(91a) 


where 


M xe i(x) 

= c^,(x) + C^ 2 (x) 

(91b) 

m x92 ( Xj 

|--c^ l (x) + c^ 2 (x) 

(91c) 

M x03 (x) 

= c“* 3 (x) - c“^ 4 (x) 

(9 Id) 

m x04 (x; 

) = c“^ 3 (x) + C^ 4 (x) 

(91e) 

( -'n|B 2 6 

B 16 ^ 12 + BJ 62 j + B 16 5 n + B 66 5 j 6 

(9 If) 

^t( ^26 

Bi/ 12 + B 66 ^ 62 j + B i6 a 16 + BggSgg 

(9 lg) 


m B,/ 12 + b 26 + b 6 / 62 


c, = 


R 


+ 2S(B 16 c 12 + B 66 5 62 + D l6 ) 
= 2 Vt 2 - S (b i6 c 12 + B 66 c 62 + D, 6 j 


(9 lh) 
(91i) 


Then, equations (4) and (7a) give 

N x0 (x) = T 0 + T,(x) - ^M xe (x) (92) 

The transverse shear stress resultant Q 0 is obtained from equation (B5); that is, 

Q b (x) = m x6 '(x) (93) 

Differentiating equations (89) and (91), and using equations (53) gives 


31 


QeW = Qho( x ) + Qbi( x ) K, + Q B 2 (x) K 2 + Q 63 (x) K 3 + Q 04 (x) K 4 


(94) 


where 


Qeo( x ) - (B| 6 *u + B 66 5 16 )[n/(x) + N x '(X)] + (B 16 « 16 + B 66 * 66 )[t/(x) + T , (x) J 

- - w/(x) _ Th 

+ (B 26 + BJ n + B 6 / 62 ) — (D | 6 + B 16 i 12 + BggCgjjWp^^x) - M x( /(x) 


Q ei (x) = - [c“vT^S^ + c“VT + s]/,(x) + [c^VT + s - c“/r - s' 


Qe 2 ( x ) = [c^vT^ - c^^tx) - [cr/TTs^ + c“/r^]/ 2 (x) 


q 63 (x) = [c^vt^ + c r.T + sy,(x) + [ c r/m - cr.T-sjf.tx) 


Q 84 ( x ) = [c“/r^ - c “vT^t?]^ 3 (x) + [c^/r^ + C r/r^]^ 4 (x) 


(95a) 

(95b) 

(95c) 

(95d) 

(95e) 


Particular Solutions 

Whenever P,(x) given by equation (29c) is nonzero, the particular solution w p (x) must be 
determined. To illustrate the solution process, consider a cylinder subjected to the surface 


tractions 

q x ( x ) = q! + q( x + q) sin(^) + ql cos(^) (96a) 

q 8 ( x ) = q“ + q« x + q; sin(^) + q B cos(^) (96b) 

q„(x) = q! + ql X + q s n sin(^) + q° cos(^) (96c) 

and the temperature-change distributions 

0 o (x) = H 0 + H[X + H s sin(^) + H c cos(^) (97a) 

©,(x) = G 0 + G, x + G s sin(^) + G c cos(^) (97b) 
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where q°, ql, qL q L q°, qi, qL qL q“, ql, qii, q», H 0 , h„ h s , h c , g 0 , g„ G s , and G c are specified 

constants. The first step in constructing the particular solution is to obtain all the input quantities 
that appear in equation (29c) for P^x). Equations (lb) and (96a) give 


NiM = -q!x-^x 2 + ^ cos(^) - 2jt 


£■*(¥) 


(98) 


and equations (7b) and (96b) give 


9eL 


T ( \ v ^ . JL o — 

i(x) = -q fl x-^-x + -jj- cos 


2jt 


sin 


(¥) 


(99) 


Equations (B17) and (97) give 


N^ h (x) = (H 0 n™ + G 0 c x h ) + (H,n™ + Gl c™)x 

+ (H s n™ + G s c x h )sin(^) + (H c n x h + G c e x h )eos(^) 


(100a) 


Likewise, 


Ne h (x) = (H o n 0 h + G 0 cD + (H in e h + G.cfjx 

+ (H s ne h + G s Ce h )sin(^) + (HpC + G c Ce”)cos(^ 


(100b) 


M™(x) = (H 0 c x h + G 0 m x h ) + (H,c™ + G,m x h )x 

+ (H,cr + G s m x h )sin(^) + (H c e™ + Gpn x h )cos(^j 


(100c) 


In addition, equations (B17), (12f), and (97) give 


~ Th , \ It t - Th _ Th \ Itt - Th _ Th \ 

T (x) = )H 0 n x6 + G 0 c xe J + (H.n^ + G[C X 0 Jx 

+ (H s n x0 h + G s c x0 h )sin(^) + (H c n™ + G c c™)eos(^) 


(lOOd) 


where 


_ Th Th JA, Th 

1^x0 — 1^x0 J^C x q 


_ Th Th U Th 

^x0 — ^x0 J^^x0 


(101a) 

(101b) 


Using equations (96c), (98), (99), and (100) with equation (29c) yields 
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(102a) 


P,(x) = Y 0 + Y t x + Y 2 x 2 

+ Y„ sin(^) + Y s2 sin(^) + Y cl cos(^) + Y c2 cos( 


where 


Yg — [q„ q x (B u 5 11 + B 16 5 16 ) q e (B n 5 16 + B 16 5 66 j 


D„ C 


+ ( n e n x (A 12 « u + A 26 # 16 ) n x6 (A l2 * 16 + A 26 * 66 )) 


d (c 0 c x (A 12 a u + A, 6 a 16 ) c xe (A 12 5 16 + A 26 « 66 )) 


Y,= 


d„ e 


qn + ( A i 2 «n + A 26 « 16 ) + p ( A 12 5 16 + A 26 a 66 j 


R' 


+ R 12*11 Y 26*16 ) n xe(Y | 2 5 16 + A 26 « 66 j ) 

(c 0 — C x (A I2 5 u + A 26 5| 6 ) — C x0 (A 12 5 16 + A 26 5 66 )) 


Gj I Th Th 

R“ 


Y,= 


l r ' 


2RD„ e 


[*1 x(Y| 2 « u + A 26 5 16 ) + q e (A l2 5 16 + A 26 5 66 )] 


Y., = 


Du e 


q n p ( n s n x (A i 2 «h + A 26 a 16 | n x0 (A 12 5 16 + A 26 « 66 )) 


+ H s ^ 2 (c x n x (B n <t n + B 16 * i6 ) n x8 (B n « l6 + B 16 * 66 )) 


+ r-. (e e c x ( A 12 a u + A 26 « 16 j c x g( A 12 *!6 + A 26 « 66 ) ] 


R 


JC / Th Th / — \ _ Th / — \ \ 

+ G s ? (rn x — c x (B u « u + B 16 « 16 j — c xe (B u 5| 6 + B 16 « 66 j J 

L/ 


Y s2 = 


( qx(Bi,«„ + B 16 « 16 ) + q 0 (B u S 16 + B 16 5 66 jj 

d„ e 

+ (qx(Ai 2 «„ + A 26 « 16 ) + q e (A 12 5 16 + A 26 « 66 J jj 


(102b) 


(102c) 


(102d) 


(102e) 


(102f) 
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(102g) 


Y, = 


l r ji 


Di! e 


[]^( ^lx(®il*ll + B 16^ 16 ) Qb(®11*16 + D 16 5 66 )) 


JlR ( < ^ x ( J ^ 12 ^ 11 ^ 26*16) "I" ^0(^12*16 + A 26*66))] 


Y,= 


D„ c 


On + ( n e n x (A 12*11 ^ 26* 16 ) 12*16 + -^26*6 


4jt 

■*■ H c 2~(^x — flx (® 11*11 "I" B 16* 16 ) — ^xe(B 11*16 + B 16 * 66 ) ) 

L (102h) 


Gc 

R 


+ ri (c e c x ( A 12 *n + A 26 « 16 ) c xS ( A 12 * 16 + A 26 « 66 ) ) 

2 

4jC / jjj yj] / — \ _ Th / — \ \ 

G c ^2 (ftR — c x [Bu*u + B 16 * 16 ) - c xB |Bn« 16 + B 16 *5j j j 


By applying the method of undetermined coefficients to equation (40) with P 0 = 0, the particular 
solution w p (x) is found to be given by 


w p (x) = W p ° + w] X + Wp X 2 

(JtX) 


+ w p sm(^J + w p sin 




— / TTX \ — 

w p cos(^) + w p cos 


(tt) 


(103) 


where 


0 



4T 


Y n - 



w„ 


4T 2 


w n 


Y^ 

4T 2 


si 




4jif 


S + 4T~ 


s2 




S + 4T 2 
L 


(104a) 

(104b) 

(104c) 


(104d) 

(104e) 
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W„ = Y 


71 

L 4 


S + 4T 2 


(104f) 


c2 




16jt 4 

L 4 



S +4T 2 


(104g) 


Examples 

Five examples are presented subsequently to illustrate the utility of the exact solution 
presented herein and some of the behavioral characteristics of laminated-composite cylinders that 
exhibit general forms of anisotropy. For each example, results are presented for length-to-radius 
ratios L/R = 0.5 and 2, and for fixed values of the radius-to-thickness ratio R/h. Cylinders with 
these two L/R values exhibit the full range of variations in behavior; that is, behavior in which 
zones of localized bending are isolated near the ends of a cylinder and behavior in which the 
bending behavior is not localized at all. Zones of localized bending are referred to herein as 
bending boundary layers. For every laminate considered, the lamina material properties are given 
by the major principal modulus E, = 18.5 x 10 6 psi, the minor principal modulus E, = 1.64 x 10 6 
psi, the shear modulus G 12 = 0.87 x 10 6 psi, the major Poisson’s ratio v 12 = 0.30, the major principal 
coefficient of thermal expansion a, = 0.25 x 10 6 /°F, and the minor principal coefficient of thermal 
expansion a 2 = 16.2 x 10 6 /°F. For all laminates considered, the ply thicknesses are identical and 
the ply thickness is taken as the total thickness h divided by the number of plies. From a practical 
standpoint, atypically thick plies are considered to be composed of groups of thinner plies with 
the same orientation. For all results presented, based on the exact solution, the Sanders-Koiter 
equations were used, which corresponds to the value u = 3/2. For several of the examples, 
simplified approximate solutions are also presented. Each simplified solution provides a 
necessary condition for establishing the correctness, and hence credibility, of the corresponding 
more complicated exact solution. Several examples also include the effects of temperature-change 
fields. For convenience, the temperature-change fields on the outside and inside of a cylinder are 

denoted by AT„ = AT(x, L) anc j Aiv = AT(x, - respectively (see equations (B 16)). 


Example 1: Self-Equilibrated End Loads 

Perhaps one of the most fundamental problems of shell theory is the attenuation behavior of 
a circular cylinder subjected to self-equilibrated end loads. For the example considered herein, the 
boundary conditions specified in the present study are 


N x (0) = 0 


(105a) 
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T(0) = 0 


(105b) 


at x = 0 and 


Qx(°) = v 0 

(105c) 

M,(0) = M„ 

(105d) 

u(L) = 0 

(106a) 

v(L) = 0 

(106b) 

w(L) = 0 

(106c) 

(3JL) = o 

(106d) 


at x = L. These boundary conditions correspond to a cylinder that is completely clamped around 
the circumference at x = L and subjected to a uniform transverse shear stress resultant and a 
uniform bending moment resultant, as depicted in figure 4. For these loading conditions, the 
particular solution given by equation (103) is identically zero-valued. 

Results were obtained for geometrically identical cylinders with two very different 

laminated-wall constructions, and for the applied loads given by M 0 = 100 lbs and V 0 = 100 lbs/ 
in. The geometries considered are characterized by the radius-to-thickness ratio R/h =100 and the 
length-to-radius ratios L/R = 0.5 and 2. In particular, results are presented in figures 4-7 for 
cylinders with laminated walls composed of all 90-degree plies and a total wall thickness h = 0.1 
in. This wall construction is referred to herein as specially orthotropic because the principal axes 
of each orthotropic ply are aligned with the x and 0 coordinate directions shown in figure 3. 
Therefore, these laminates exhibit no form of anisotropy associated with coupling between 
extension and shear, coupling between bending and twisting, or coupling between any membrane 
deformation and any bending deformation. As a result, the cylinders with this wall construction 
do not exhibit any circumferential displacements. The values of the nonzero constitutive matrices 
computed for this laminate are given by 


A ii A 12 A i 6 


AAA 

JT. 12 JT- 22 -^26 

- 

A 16 A 26 A 66 



0.165319 

0.0495957 

0 


0.0495957 0 

1.86488 0 

0 0.087000 


io 6 !^ 


in. 


(107a) 
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Dll Di: D,6 


0.137766 

0.0413297 

0 

3 


D 12 D 2 2 0,6 

= 

0.0413297 

1.55407 

0 

x 10 in. -lbs 

(107b) 

Dl 6 D 26 D 6 6 


0 

0 

0.0725000 




Similarly, the constitutive expressions defined by equations (15) 






*u a l6 


6.04891 

0 

*16 *66 


0 

1 1 .494 


x 10 


in. 

lbs 


(108a) 


/M 

/ 


1 0.3 1 
\ 0 / 


(108b) 


(108c) 


The values of S and T, defined by equations (27) and (39), respectively, that are used to 
determine the basis functions of the solution space are given by S = 0 and T = 5.79409 in.' 2 . 

In contrast to the specially orthotropic wall, results are presented in figures 8-13 for cylinders 
with laminated walls composed of all 10-degree plies and a total wall thickness h = 0.1 in. This 
wall construction exhibits anisotropy in the form of coupling between membrane expansion or 
contraction and shearing deformations, and coupling between pure bending and twisting 
deformations. The values of the nonzero constitutive matrices computed are given by 


A i] A 12 A 16 

AAA 

12 -^22 26 

A 16 A 26 A 66 


1.76734 

0.0958898 

0.272513 


0.0958898 0.272513 

0.170273 0.0181288 

0.0181288 0.133294 


x io 6 !^ 
in. 


(109a) 


D„ D 12 D 16 


D 12 D 22 D 26 

= 

D 16 D26D66 



1.47278 

0.0799082 

0.227094 


0.0799082 0.227094 

0.141894 0.0151073 

0.0151073 0.111078 


3 

x 10 m.-lbs 


(109b) 


Similarly, the constitutive expressions defined by equations (15) are given by 

*11 a l6 
*16 *66 


0.826304 - 1.6893 
- 1.6893 10.9577 


x 10 


in. 

lbs 


(110a) 



/ 0.0486091 \ 
\ 0.0366266 / 


(110b) 
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(110c) 


/ «i2 \ _ / 0.575447 \ 
\* 62 f~ \ 3.731960 / 


x 10 4 in. 


For this anisotropic-wall construction, S = 4.23574 x 10 5 in.' 2 and T = 0.529147 in.' 2 . 

M 

Graphs of the normalized radial displacement w/h, bending stress resultant -~p transverse- 

Mo 

O N 

shear stress resultant 7K and circumferential stress resultant for the cylinders with all 90- 

V„ v 0 

degree plies are shown in figures 4-7, respectively, as a function of the axial distance x normalized 
by the cylinder length L. Two curves appear in each figure that correspond to values of L/R = 0.5 
and 2. These results exhibit the common trend in which the effects of the edge loads attenuate as 
x/L increases, with the longer cylinder exhibiting the larger attenuation rate. Likewise, graphs of 
the normalized axial displacement u/h, circumferential displacement v/h, radial displacement 

w/h, bending stress resultant transverse-shear stress resultant 8\ and circumferential 

M 0 V 0 

stress resultant ^' !l for the cylinders with all 10-degree plies are shown in figures 8-13, 

Vo 

respectively, as a function of the normalized axial distance x/L. The results for these cylinders 
with L/R = 2 also exhibit the expected response-attenuation behavior. It is noteworthy to point out 
that the response quantities for the shorter cylinder with all 10-degree plies do not attenuate as x/ 
L increases, unlike the corresponding cylinders with all 90-degree plies. This difference in 
attenuation behavior is due to the higher axial stiffness of the cylinders with all 10-degree plies. 
The results in figure 9 show relatively small circumferential displacements, as expected, that are 
generated by the wall anisotropy. For both wall constructions, circumferential stresses (figure 13) 
are generated that are significantly larger in magnitude than the corresponding value of the 
applied transverse-shear load. Moreover, the largest magnitude of the circumferential stresses in 
the cylinders with all 90-degree plies (figure 7) are approximately six times larger than those in 
the cylinders with all 10-degree plies. 


Example 2: Uniform Axial End Load 

Another fundamental problem of shell theory is the response-attenuation behavior of a 
circular cylinder subjected to a uniform axial end load. This problem is particularly important in 
the study of cylinder buckling. For the corresponding example considered in the present study, 
the boundary conditions are given by 


u(0) = 0 

(Hla) 

vfO) = 0 

(111b) 
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w(0) = 0 


(lllc) 


M x (0) = 0 

(11 Id) 

N,(L) = N l 

(112a) 

V (L) = 0 

(112b) 

w(L) = 0 

(112c) 

M X (L) = 0 

(112d) 


at x = L. For these loading conditions, the particular solution given by equation (103) is 
identically zero-valued. 

Results for this example were also obtained for cylinders with radius-to-thickness ratio R/h 
= 100 and the length-to-radius ratios L/R = 0.5 and 2. Specifically, results are presented in figures 
14-19 for cylinders with laminated walls composed of two inner 90-degree plies, two outer 10- 
degree plies, and a total wall thickness h = 0. 1 in. This wall construction exhibits the full extent 
of anisotropy possible in the exact solution; that is, coupling between membrane expansion or 
contraction and shearing deformations, coupling between pure bending and twisting 
deformations, and coupling between membrane and bending-twisting deformations. The values 
of the constitutive matrices that were computed are given by 


An A 12 A i 


16 A 26 A 66 


0.966328 

0.0727428 

0.136257 


0.0727428 

1.01758 

0.00906441 


0.136257 

0.00906441 

0.110147 


xlO 6 ^ 

in. 


(113a) 


B , , 

B |2 B | 6 


2.00252 

0.0578677 

0.340641 

B, 2 

B 22 B 26 

= 

0.0578677 

- 2.11826 

0.022661 

B | 6 

B 26 B 66 


0.340641 

0.022661 

0.0578677 


x 10 lbs 


(113b) 


D„ D 12 D ,6 


Dl 2 D 2 2 D 26 

= 

Di 6 D 26 D 66 



8.05273 

0.60619 

1.13547 


0.60619 

8.4798 

0.0755367 


1.13547 

0.0755367 

0.917892 


2 

x 10 in. -lbs 


(113c) 


Similarly, the constitutive expressions defined by equations (15) are given by 
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(114a) 


< 2-11 < Z l6 


1.25506 

- 1.5559 

^16 ^66 


- 1.5559 

10.9932 


x 10 


in. 

lbs 



/ -0.0771401 \ 
\ 0.0131601 / 


/ c i2 \ 
\«62 / 


/ 1.98063 \ 
\ 0.64773 / 


x 10 " in. 


For this laminate, S = 0.119157 in.' 2 and T = 2.55869 in.' 2 . 


(114b) 


(114c) 


Graphs of the normalized axial displacement u/h, radial displacement w/h, bending stress 

resultants anc j M 0 twisting stress resultant M - e and circumferential stress resultant ^ 
N L h N L h’ N L h’ N l 

are shown in figures 14-19, respectively, as a function of the axial distance x normalized by the 
cylinder length L. Two pairs of blue and red curves appear in each figure that correspond to values 
of L/R = 0.5 and 2. The red curves correspond to an approximate solution, given in Appendix C, 
that is based on the equations of linear membrane theory, in which the bending boundary layer is 
ignored such that the radial displacement is presumed to be constant valued. In the figures with 
only one red curve, the curve corresponds to both values of L/R. For all results, a value of 

N l = 1000 lbs/in. was used for the end load. The blue curves correspond to exact solution 
presented herein. 

The results in figure 14 show a nearly linear variation in the axial displacement, with the 
largest displacements being exhibited by the longer cylinder. Moreover, the results indicate that 
the approximate solution yields nearly the same results as the exact solution for both values of 
L/R. The normalized axial displacements at x = L were computed as u/h = 0.216 and 0.060 for 
L/R = 2 and 0.5, respectively. In contrast, the results in figures 15-19 show the presence of 
bending boundary layers in the longer cylinders with L/R = 2 that attenuate to the approximate 
solution, as expected, at distances from the cylinder ends that are about 10% of the cylinder 
length. The results also indicate that most of the the response quantities for the shorter cylinders 
with L/R = 0.5 do not exhibit zones of bending attenuation and, as a result, the approximate 
solution is inadequate for estimating the interior response of the cylinders. The normalized radial 
displacements at x = L/2 were computed as w/h = -0.007 and -0.004 for L/R = 2 and 0.5, 
respectively. 


Example 3: Uniform Temperature and Internal Pressure Fields 


A very practical problem of shell theory deals with the behavior of a circular cylinder 
subjected to uniform internal pressure and different, uniform internal and external temperature 
fields. For the corresponding example considered herein, the boundary conditions are given by 
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u(0) = 0 


(115a) 


at x = 0 and 


v(0) = 0 

(115b) 

w(0) = 0 

(115c) 

|3°(°) = 0 

(115d) 

u(L) = 0 

(116a) 

v(L) = 0 

(116b) 

w(L) = 0 

(116c) 

P°(L) = 0 

(116d) 


at x = L. Thus, the cylinder is completely clamped around the circumference at each end. The 
applied internal pressure given by q„M = q" = too psi and the applied thermal load is given by an 
external temperature increase AT 0 = 300 "F and an internal temperature increase AT = 100 °F. 
From equations (B16) and (97) it follows that 0 o (x) = H 0 = 200 °F and 0,(x) = G 0 = 2000 °F/in. 
For these loading conditions, the particular solution given by equation (103) reduces to 


— , , — ° 
w p (x) = W = J 

4T 


(117a) 


where Y 0 given by equation (102b) reduces to 


Yn = 


l r o , H 


d„ e 


[fin ■*" ( n e n x ( A 12 « n + A 26 a 16 ) n xB (A 12 5 16 + A 26 « 66 )) 


+ 


G 


( c e c x (A 12 a n + A 26 5 16 ) c x 6 (A 12 a 16 + A 26 5 66 j j 


(117b) 


For this example, results were obtained for cylinders with radius-to-thickness ratio R/h = 100 
and the length-to-radius ratios L/R = 0.5 and 2. Selected results are presented in figures 20-27 for 
the previous generally anisotropic cylinders with laminated walls composed of two inner 90- 
degree plies, two outer 10-degree plies, and a total wall thickness h = 0.1 in. The values of the 
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constitutive matrices computed are given by equations (113) and (114) and, as before, S = 

0. 1 19 157 in.' 2 and T = 2.55869 in.' 2 . The values of the thermal constitutive parameters in equation 
(B17) that were computed are given by 


| 2 S7 ) j»- 

(-0.121494) F-in. 


(118a) 


( - 0.0172256 ) ,, 
0.0172256 } ^ 
( -0.00303735 J F 


(118b) 


m. 


m n 


m. 


H 


0.00166795 \ 

0.00163225 
0.000101245 j 


in. -lbs 


(118c) 


Graphs of the normalized axial displacement u/h, circumferential displacement v/h, radial 

N 

displacement w/h, circumferential stress resultant — — - — bending stress resultants 

q„R + N c 

, „ M * Th , and 7 „ Mb Th . , twisting stress resultant , , Mxe T . , , and transverse-shear 
( q.R + N™)h (q':R + N:)h’ 8 (qS + NDh’ 

stress resultant w are shown in figures 20-27, respectively, as a function of the axial 

q„R + N 0 

distance x normalized by the cylinder length L. Two pairs of blue curves and coincident red lines 
appear in each figure that correspond to values of L/R = 0.5 and 2. The red lines correspond to an 
approximate solution given in Appendix D that is based on the equations of linear membrane 
theory, in which the bending boundary layer is ignored such that the radial displacement is 
presumed to be constant valued. 

The results in figures 20 and 21 show relatively small-magnitude antisymmetric distributions 
of the axial and circumferential displacements that attenuate to a value of zero in the central region 
of the cylinders, with the larger magnitudes occuring near the ends. The corresponding 
displacements predicted by the approximate solution are zero valued, as indicated by the red line 
in each figure. The exact solution gives the constant values N x = -288.68 lbs/in. and -270.10 lbs/ 
in. for the cylinders with L/R = 0.5 and 2, respectively. In addition, the exact solution gives the 

constant values T = 40.15 lbs/in. and 42.47 lbs/in. for the cylinders with L/R = 0.5 and 2, 

respectively. The approximate solution gives N x = -263.90 lbs/in. and f =43.25 lbs/in., 
regardless of the value for L/R. 
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The results in figures 22-27 show the presence of bending boundary layers in the longer 
cylinders with L/R = 2 that attenuate to the approximate solution, as expected, at distances from 
the cylinder ends that are about 10% of the cylinder length. The results also indicate that most of 
the response quantities for the shorter cylinder with L/R = 0.5 exhibit small centrally located 
zones of attenuation and, as a result, the approximate solution is inadequate for estimating the 
interior response of these cylinders. The approximate solution yields the values w/h = 0.14, 


N e 

q!R + Nl 


= 0.701, 


M v 


(q^R + Ng h )h 


= 0.006, 


M„ 


(qnR + N™ )h 


= - 0.255, and 


Th 

7 — xfl T , , = 0.008 with N e =426.19 lbs/in. The values agree well with the exact solution 

(q!R + N B )h 

in the central, attenuated region of the longer cylinders. However, even in the central region of the 
shorter cylinders, there are significant differences between the approximate and exact solutions 
for several of these quantities. 


Example 4: Linearly Varying Temperature and Internal Pressure Fields 

This example corresponds to cylinder proportions and loading conditions that are typical of 
a large-scale launch vehicle. Specifically, a circular cylinder subjected to an axial end load, an end 
circumferential shear load, a linearly varying internal pressure field, associated with a liquid fuel, 
and different, linearly varying internal and external temperature fields. For this example, the 
boundary conditions specified in the present study are 


at x = 0 and 


N x (0) = N 0 

(119a) 

f (0) = T 0 

(119b) 

w(0) = 0 

(119c) 

M x (0) = 0 

(119d) 

u(L) = 0 

(120a) 

v(L) = 0 

(120b) 

w(L) = 0 

(120c) 
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M x (L) = 0 


(120d) 


at x = L. The values of the applied end loads used in the present example are given by N„ = 1000 

lbs/in. and T„ = 200 lbs/in. The applied pressure load is given by q„M = p„ + y(L - x), where p u 
= 30 psi represents an ullage pressure and y = 0.04123 lb/in. 3 is the specific weight of the fluid 
inside the cylinder. Thus, equation (96c) reduces to q„M = q“ + q! x, where q“ = p u + yL and 

q! = - y. The applied thermal load is given by an external temperature rise AT 0 (x) = - 500 + 100 

x/L °F and an internal temperature rise AT.(x) = - 500 + 50 x/L °F. From equations (B 16) and (97) 

it follows that, 0 o (x) = H o + H,x and 0 1 (x) = G o + G 1 x, where H 0 = - 500 °F, G 0 = 0 “F/in., 

H [ = 75/L °F/in., and G, = 50/(hL) °F/in. 2 . For these loading conditions, the particular solution 

_ o i _o Y 0 _i Y, 

given by equation (103) reduces to W P ( X ) = w p + w p x, where W P = pp and w p = — i . in 

addition, Y 0 and Y p given by equations (102), reduce to 


Y = 


l r “ , H 0 


d„ e 


[fin + ( n e n x (A i 2 '* i i + A 26 5 16 ) n x0 ( A 12 a 16 + A 26 a 66 ) ) 


Go 

R 


(^e c x (A 12 *n + A 26 # 16 j c x e( A 12 a 16 + A 26 <* 66 j | 


(121a) 


and 


A i — [q n + p (n e n x (A 12 5 u + A 26 5 16 ) n xB (A| 2 « 16 + A 26 S 66 jj 


d u e 


R 


"t" (^0 c x ( A 12 «t u + A 26 <* 16 j c xB (A 12 « 16 + A 26 <* 66 ) ) 


(121b) 


For this example, results were obtained for cylinders with radius-to-thickness ratio R/h = 330 
and the length-to-radius ratios L/R = 0.5 and 2. Selected results are presented in figures 28-35 for 
the previously considered generally anisotropic cylinders with laminated walls composed of two 
inner 90-degree plies, two outer 10-degree plies, but with a total wall thickness h = 0.5 in. The 
values of the constitutive matrices computed are given by 


A a A in A 16 


AAA 

J-V 12 -^22 xx 26 

= 

A 16 A 26 A 66 



4.83164 

0.363714 

0.681283 


0.363714 

5.08788 

0.045322 


0.681283 

0.045322 

0.550735 


x 10 


lbs 


in. 


(122a) 
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B n B 12 B i6 


5.00630 

0.144669 

0.851603 

B 12 B 22 B 26 

= 

0.144669 

- 5.29564 

0.0566525 

B| 6 B 26 B 66 


0.851603 

0.0566525 

0.144669 


(122b) 


D n 

D,2 

D 16 


100659 

7577.37 

14193.4 

Dij 

d 22 

D 2 6 

= 

7577.37 

105997 

944.209 

Di« 

C>26 

D “. 


14193.4 

944.209 

11473.7 


in. -lbs 


(122c) 


and the constitutive expressions defined by equations (15) are given by 


#11 # j 6 


2.50793 

- 3.10445 

^16 ^66 


- 3.10445 

2.19917 


x 10 


in. 

lbs 


(123a) 



/ -0.0771309 \ 
\ 0.0131291 / 


(123b) 


/ 5 12 \ = _ / 0.0990770 \ : 
\ e 62 j \ 0.0321475 / 


(123c) 


The values of S and T, defined by equations (27) and (39), respectively, that are used to 
determine the basis functions of the solution space were found to be given by S = 0.001444 in.' 2 
and T = 0.031012 in.' 2 , and the values of the thermal constitutive parameters in equation (B 17) 
are given by 


I 10.0077 ) ,, 

{ 9.79348 } 

( - 0.607469 j F- 


in. 


(124a) 


- 0.430641 
0.430641 
-0.0759336 


lbs 


(124b) 


m. 


m 


m. 


I 0.208494 ) . ,, 

| 0.204031 ) in - 0 lbs 

\- 0.0126556 ) F 


(124c) 


Graphs of the normalized axial displacement u/h, circumferential displacement v/h, radial 
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displacement w/h, circumferential stress resultant 9 . , bending stress resultants 

(p u + yL)R 


t - r — and — , twisting stress resultant — , and transverse-shear 

(p u + Y L)Rh (p u + yL)Rh 6 ( Pll + yL)Rh 

stress resultant are shown in figures 28-35, respectively, as a function of the axial 

(Pu + yL)R 

distance x normalized by the cylinder length L. The denominator term p u + yL represents the 
maximum pressure in the cylinder. Two pairs of blue and red curves appear in each figure that 
correspond to values of L/R = 0.5 and 2. The red curves correspond to an approximate solution 
given in Appendix E that is found by applying the equations of linear membrane theory, ignoring 
the bending boundary layer and presuming that the radial displacements vary like the pressure and 
temperature fields. In the figures with only one red curve, the curve corresponds to both values of 
L/R. 


The results in figures 28 and 29 show a nearly linear variation in the axial and circumferential 
displacements, respectively, with the largest displacement magnitudes being exhibited by the 
longer cylinder. Moreover, the results indicate that the approximate solution yields nearly the 
same results as the exact solution for both values of L/R. In contrast, the results in figures 30-35 
show the presence of bending boundary layers that attenuate to the linearly varying approximate 
solution at distances from the cylinder ends that are about 10% and 20% of the cylinder length for 
L/R = 2 and 0.5, respectively. The linear variation in the response quantities is generally more 
pronounced in the longer cylinders because the pressure increases with increasing cylinder length. 


Example 5: Sinusoidal Pressure and Temperature Fields 

This example illustrates the behavior of a circular cylinder subjected to sinusoidal internal 
pressure and different, sinusoidal internal and external temperature fields. The cylinder is fixed at 
one end and free to expand axially and twist circumferentially at the other end. The corresponding 
boundary conditions specified in the present study are 


u(0) = 0 

(125a) 

v(0) = 0 

(125b) 

w(0) = 0 

(125c) 

p°(0) = 0 

(125d) 
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at x = 0 and 


N x (0) = 0 

(126a) 

o 

ll 

o' 

<H 

(126b) 

w(L) = 0 

(126c) 

P°(L) = 0 

(126d) 


at x = L. The applied mechanical load given by q„M = q“ sin(^), where q! = 1000 psi. The 
applied thermal load is given by an external temperature rise AT 0 (x) = 300 sin(jtx/L) °F and an 
internal temperature rise AT(x) = 100 sin(jtx/L) °L. Equations (B16) and (97) give 0 o (x) = H 
sin(jtx/L) and 0,(x) = G s sin(itx/L) where H s = 200 °F and G s = 200/h "F/in. For these loading 
conditions, the particular solution given by equation (103) reduces to 

w p (x) = w;‘ sin(™) (127) 


where w' 1 is given by equation (104d). 

For this example, results were obtained for cylinders with radius-to-thickness ratio R/h = 100 
and the length-to-radius ratios L/R = 0.5 and 2. Selected results are presented in figures 36-45 for 
generally anisotropic cylinders with laminated walls composed of two inner 90-degree plies, two 
outer 10-degree plies, and a total wall thickness h = 0. 1 in. The values of the constitutive matrices 
are given by equations (113) and (114), and the values of the thermal constitutive parameters in 
equation (B17) are given by equations (118). Like before, S = 0.119157 in. 2 andT = 2.55869 in.' 2 . 

Graphs of the normalized axial displacement u/h, circumferential displacement v/h, and 
radial displacement w/h are shown in figures 36-38, respectively, as a function of the axial 

distance x normalized by the cylinder length L, and the axial rotation (h is shown in figure 39 as 

a function of x/L. Graphs of the normalized circumferential stress resultant M, bending stress 

N 0 

resultants Mi anc j Mi twisting stress resultant Mi and transverse-shear stress resultants 
M x M e ’ M x 

MM anc j MM are shown in figures 40-45, respectively, as a function of x/L. In these 
M x M x 

nondimensional quantities, N e is given by equation (LlOb), M x is given by equation (L7b), and 
M„ is given by equation (LI lb). One blue curve and one red curve is shown in each part of the 
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two-part figures. The red curves correspond to a simplified exact trigonometric solution, given in 
Appendix F, for the corresponding cylinder with simply supported ends instead of clamped ends. 
The blue curves represent the exact solution obtained herein. 

The results in figures 36-45 show that the difference between the exact solution (blue 
response curves) for the cylinders with rotation restrained at the ends and the exact trigonometric 
solution for simply supported ends is much smaller for the cylinders with L/R = 2 than for those 
with L/R = 0.5, as expected. For the longer cylinders, significant differences between most of the 
corresponding response curves appear to be confined to edge zones that are about 10% of the 
cylinder length. For the shorter cylinders, the zones of significant difference grow to about 20% 
of the cylinder length. For both cylinder lengths, the greatest overall discrepancies between the 
exact and corresponding trigonometric solutions are exhibted by the axial displacements and the 
bending stress resultant M x . 


Concluding Remarks 

A detailed exact solution has been presented for generally laminated circular cylinders that 
undergo axisymmetric deformations. The cylinders exhibit the full range of anisotropy that is 
possible in first-approximation shell theories that neglect high-order effects such as transverse- 
shear deformation. The overall solution has been formulated in a general, systematic way and is 
based on the solution of a single fourth-order, nonhomogeneous ordinary differential equation 
with constant coefficients in which the radial displacement is the dependent variable. Moreover, 
positive-definiteness of the strain energy has been used to define uniquely the form of the basis 
functions spanning the solution space of the ordinary differential equation. This approach appears 
to be a unique feature of the analysis presented herein. 

Axisymmetric loading conditions have been considered that include edge loads, surface 
tractions, and temperature fields. In addition, all possible corresponding boundary conditions 
have been considered, and several particular solutions to the fourth-order, nonhomogeneous 
ordinary differential equation governing the response have been presented. Results have been 
presented for five examples that demonstrate a wide range of behavior. Several of the results 
illustrate how laminated-wall construction and cylinder length affect the size of localized zones 
of bending deformations that occur near the cylinder ends. For some of the examples, the results 
indicate highly localized bending deformations that attenuate rapidly. For other examples, the 
bending deformations exhibit very little attenuation. In addition, approximate solutions for 
several of the examples have been presented that add credibility to the corresponding results 
obtained from the exact solution and that provide insight into the possible pitfalls of making 
certain types of simplifications in a preliminary-design process. For example, several results are 
presented that illustrate conditions for which sizing a laminated-composite cylinder based on an 
approximate linear membrane analysis is inappropriate. 

Results obtained from the exact solution presented herein also have the potential to serve as 
fundamental benchmarks for the validation of numerical solution methods such as the finite- 
element method. Likewise, results obtained from the exact solution have the potential to provide 


49 


insight into the extent of discretization that might be needed for analyzing accurately shells with 
zones of localized bending caused by cutouts, stiffeners, or other types of discontinuities and 
gradients. 
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Figure 1. Space Shuttle external tank components. 
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Figure 2. Conceptual lunar lander. 



Figure 3. Cylinder geometry, surface coordinate system, and lamina fiber orientation. 


54 



Figure 4. Radial displacement of laminated cylinder with all 90-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 



Figure 5. Bending stress resultant for laminated cylinder with all 90-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 
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Figure 6. Transverse shear stress resultant for laminated cylinder with all 90-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 



Figure 7. Circumferential stress resultant for laminated cylinder with all 90-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 
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Figure 8. Axial displacement of laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 



Figure 9. Circumferential displacement of laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 
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Figure 10. Radial displacement of laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 



Figure 11. Bending stress resultant for laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 
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Figure 12. Transverse shear stress resultant for laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 



Figure 13. Circumferential stress resultant for laminated cylinder with all 10-deg plies, subjected to edge loads at x = 0, 
and clamped at x = L (R/h = 100). 
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Figure 14. Axial displacement of [90 2 /10 2 ] T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 15. Radial displacement of [90 2 /10 2 ] T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 16. Bending stress resultants of [90 2 /10J T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 17. Bending stress resultants of [90 2 /10,] T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 18. Twisting stress resultants of [90 2 /10 2 ] T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 19. Circumferential stress resultants of [90 2 /10 2 ] T laminated cylinder subjected to an axial edge load at x = L, 
with u, v, and w fixed at x = 0, and with v and w fixed at x = L (R/h = 100). 
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Figure 20. Axial displacement of [90 2 /10 2 ] T laminated cylinder clamped at both ends and subjected internal 
pressure and uniform interior and exterior temperature fields (R/h = 100). 



Figure 21. Circumferential displacement of [90 2 /10 2 ] T laminated cylinder clamped at both ends and subjected internal 
pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 22. Radial displacement of [90 j/10J t laminated cylinder clamped at both ends and subjected internal 
pressure and uniform interior and exterior temperature fields (R/h = 100). 



Figure 23. Circumferential stress resultant of [90 2 /10J T laminated cylinder clamped at both ends and 
subjected internal pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 24. Bending stress resultant of [90 2 /1 0 2 ] , laminated cylinder clamped at both ends and subjected to 
internal pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 25. Bending stress resultant of [90 2 /10 2 ] T laminated cylinder clamped at both ends and subjected to 
internal pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 26. Twisting stress resultant of [90 2 /10J T laminated cylinder clamped at both ends and subjected to 
internal pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 27. Transverse-shear stress resultant of [90 2 /10 2 ] T laminated cylinder clamped at both ends and 
subjected to internal pressure and uniform interior and exterior temperature fields (R/h = 100). 
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Figure 28. Axial displacement of [90 2 /10 2 ] T laminated cylinder subjected to edge loads at x = 0, linearly varying 
pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 29. Circumferential displacement of [90 2 /10J T laminated cylinder subjected to edge loads at x = 0, linearly 
varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 30. Radial displacement of [90 2 /10J T laminated cylinder subjected to edge loads at x = 0, linearly 
varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 31. Circumferential stress resultant of [90 2 /10 2 ] T laminated cylinder subjected to edge loads at x = 0, 
linearly varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 32. Bending stress resultant of [90 2 /10J T laminated cylinder subjected to edge loads at x = 0, linearly 
varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 33. Bending stress resultant of [90 2 /10 2 ] T laminated cylinder subjected to edge loads at x = 0, linearly 
varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 34. Twisting stress resultant of [90 2 /10 2 ] T laminated cylinder subjected to edge loads at x = 0, linearly 
varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 



Figure 35. Transverse-shear stress resultant of [90 2 /10,] T laminated cylinder subjected to edge loads at x = 0, 
linearly varying pressure and temperatures, with u, v, and w fixed at x = L, and with w fixed at x = 0 (R/h = 330). 
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Figure 36. Axial displacement of [90 2 /10 2 ] T laminated cylinder subjected to sinusoidal pressure and temperatures, 
fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential motion (R/h = 100). 


71 






Axial distance, x/L 

(b) Length-to-radius ratio, L/R = 2 



(a) Length-to-radius ratio, L/R = 0.5 


Figure 37. Circumferential displacement of [90,/10 2 ] T laminated cylinder subjected to sinusoidal pressure 
and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = fOO). 
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Figure 38. Radial displacement of [90 2 /10 2 ] T laminated cylinder subjected to sinusoidal pressure 

and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 

motion (R/h = 100). 
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Figure 39. Axial rotation of [90 2 /10 2 ] T laminated cylinder subjected to sinusoidal pressure 

and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 

motion (R/h = fOO). 
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Figure 40. Circumferential stress resultant of [90 2 /10 2 ] T laminated cylinder subjected to sinusoidal pressure 
and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 


75 






15 


Bending stress 
resultant, 

M* 

M x 


10 



.5 L — . . — ■ — ■ . — * . 

0 0.2 0.4 0.6 0.8 1.0 

Axial distance, x/L 

(b) Length-to-radius ratio, L/R = 2 


Bending stress 
resultant, 

Mx 

M x 



(a) Length-to-radius ratio, L/R = 0.5 


Figure 41. Bending stress resultant of [90 2 /10 2 ] T laminated cylinder subjected to sinusoidal pressure and 
temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 
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Figure 42. Bending stress resultant of [90 2 /10J T laminated cylinder subjected to sinusoidal pressure and 
temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 
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Figure 43. Twisting stress resultant of [90 2 /10J T laminated cylinder subjected to sinusoidal pressure and 
temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 
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Figure 44. NTransverse-shear stress resultant of [90,/10 2 ] T laminated cylinder subjected to sinusoidal pressure 
and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 
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Figure 45. Transverse-shear stress resultant of [90 2 /10,] T laminated cylinder subjected to sinusoidal pressure 
and temperatures, fixed at x = 0, and partially fixed at x = L to permit unrestrained axial and circumferential 
motion (R/h = 100). 
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Appendix A 

Sanders-Koiter Equations 

The linear Sanders-Koiter shell equations 24,25 are presented in this appendix for a right- 
circular cylinder with a radius that is given by R. For these equations, x is the axial coordinate, 0 
is the circumferential coordinate, and z is the radial coordinate, as shown in figure 3. First, the 
equilibrium equations are presented, then the kinematic equations and the constitutive equations 
are presented. Last, the boundary conditions are given for a complete right-circular cylinder at an 
edge that is given by a constant value of x. 


Equilibrium Equations 

The equilibrium equations are given in a form similar to those found in Ref. 24; that is, 


3N X i 0N xe c 2 0M X0 

0x R 00 9R 2 00 4x 

(Al) 

dN x0 l 0N B c ln c 2 dM x0 

dx R 00 R^ 9 2R 0x + 49_u 

(A2) 

0Q* 1 0Qo N e + = o 

0x R 00 R 4 " 

(A3) 

0M X i 0M xa 
0x R 00 ^ x 

(A4) 

0M xB i 0M O 
0x R 00 ^ 0 

(A5) 


where N x , N fl , and N x6 are the membrane stress resultants; Q x and Q 0 are the transverse shear- 
stress resultants; M x , M 0 , and M x0 are the bending stress resultants; q x , q 0 , and q n are the applied 
surface tractions; and c, and c, are constants that identify the equations of other shell theories that 
are considered herein. In particular, the Sanders-Koiter equations are given by c, = c 2 = 1 and the 
Love-Kirchhoff equations are give by c, = 1 and c 2 = 0. Donnell’s equations are given by c, = c, 
= 0. This convention is used throughout the present study. 


Kinematic Equations 


The kinematic equations are given by 

8 " 

°X 


du 

5x 


(A6) 
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(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 


o _ i ap° e _ c, av i a 2 w 
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(A13) 
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2 a 2 w u j_ \ av _ c i au 

Raxae Rl 1 2^ 2 jax 2 r 230 


(A 14) 


where u, v, and w are the axial, circumferential, and radial displacements of a point of the shell 
middle surface; £° , £ 0 , and Y°e are the membrane strains; P° , P° e , and P“ are the rotations; and k° , 

K° e , and K° e are the bending strains. The displacement w is positive when it is outward from the 
cylinder reference surface. 


Constitutive Equations 

The isothermal constitutive equations are given in matrix form by 


( N x \ 
N 0 
) N x0 

- 

A a A 12 A 16 

A 12 A 22 A 26 
A 16 A 26 A 66 

B j, B 12 B j 6 
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D.l D 12 D i6 
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19 12 D 22 1^26 
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k b 

\ M x0 j 


b 16 b 26 b 66 

D,6 d 26 d 66 


O 

y ^x8 


( Th \ 

1 N x 


N 
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M x 
M e 

, it 71 

V M x0 J 


(A15) 
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where the subscripted A, B, and D terms of the matrix are the stiffnesses of laminated composite 
shells that are obtained from the Love-Kirchhoff shell theory. Moreover, the constitutive terms 
in Eq. (A 15) are identical to those for laminated-composite plates that are given in Ref. 34 on pp. 
198 and 243-244. The terms with the superscript "Th" are the ficticious thermal stress resultants. 


Boundary Conditions 

The boundary conditions for an edge that is defined by a constant value of the axial 
coordinate x are given by 


N x = N(0) or u = u(0) 


(A16) 


N + — 

1 ^x0 “ 



|M xe = T(0) or v = v(0) 


(A 17) 



(A 18) 


M x = M(0) or (3° = (3(0) 


(A 19) 


where u(0), v(0), and w(0) are applied edge displacements; (3(0) is an applied edge rotation; and 
N(0), T(0), V(0), and M(0) are applied edge loads. 
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Appendix B 

Equations for Axisymmetry 


The linear Sanders-Koiter shell equations that are presented in Appendix A for a right- 
circular cylinder with a radius R are specialized in this appendix for the case of axisymmetric 
behavior. For these equations, x and 0 denote the axial and circumferential coordinates, 
respectively. The specialization to axial symmetry is conducted by eliminating all terms in the 
equations of Appendix A that are differentiatied with respect to the circumferential coordinate, 0. 
First, the equilibrium equations, the kinematic equations, and the constitutive equations are 
presented. Then, the boundary conditions are given for a complete right-circular cylinder at an 
edge that is given by a constant value of x. 


Equilibrium Equations 

The equilibrium equations for axisymmetric behavior are given by 


dN x . 

dx + 0-W = ° 

(Bl) 

c,~ c 2 dM x0 , . n 

+ R Qh+ 2R dx + 4e(x)-0 

(B2) 

dQ x N e , . _ 

dx - R + «■<*>= 0 

(B3) 

d ^ Ix Q x = 0 

dx 

(B4) 

d ^ IxB Q« = 0 

dx 

(B5) 


where the membrane stress resultants N , N„, and N ■ the transverse shear-stress resultants Q and 

x’ 0’ X0’ ^-x 

Q e ; the bending stress resultants M x , M 0 , and M x0 ; and the applied surface tractions q x , q fl , and q n 
are functions of only the axial coordinate, x. 


Kinematic Equations 

The kinematic equations are given by 


p o du 

dx 



(B6) 

(B7) 
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(B8) 

p: = - 4^ 
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(B 11) 

dp x d w 

X_ dx“ dx ^ 

(B12) 

K° 0 = O 

(B13) 

1 B° + dJ3e - 1 (c + ! c ) dv 
R 1 " dx Rl c, + 2 c ^dx 

(B14) 


where the middle-surface displacements u, v, and w; the membrane strains e°, eL and Y° e ; the 

rotations p", P° 0 , and p°; and the bending strains k", k° b , and k ° 0 are functions of only the axial 
coordinate, x. 


Constitutive Equations 

The constitutive equations reduce to 
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For thin shells, it is convenient to introduce a temperature-change function 


AT(x, z) = 0„(x) + z 0,(x) 


such that 


(B15) 


(B16a) 
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(B 16b) 


and 


0 o (x) = i[AT(x, |) + AT(x, - |)] 
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(B17) 


where the two vectors on the right-hand- side of the equal sign are the constants that result from 
substituting equation (B16) into expressions for the ficticious thermal stress resultants, like those 
given on pp. 243-244 of reference 34. 


Boundary Conditions 

The boundary conditions for an edge that is defined by a constant value of the axial 
coordinate, x, are given by 


N x = N 

or u = u 

(B18) 

|c, + 2^2 

i> 

ll 

> 

S-H 
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ll 

® 

(B19) 

Qx = v 

or w = w 

(B20) 

M x = M 

or P° x = (3 

(B21) 


where the applied edge displacements u, v, and w; the applied edge rotation P; and the applied 
edge loads N, T, V, and M are all constants. 
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Appendix C 

Approximate Solution for Uniform Axial End Load 

For this example, the boundary conditions are specified as 


at x = 0 and 


u(0) = 0 (Cla) 

v(o) = 0 (Clb) 

w(0) = 0 (Clc) 

M x (0) = 0 (Cld) 

n x (l) = n l (C2a) 

v(L) = 0 (C2b) 

w(L) = 0 (C2c) 

M X (L) = 0 (C2d) 


at x = L. The only applied load is that defined by equation (C2a). An approximate solution is 
found by applying the equations of linear membrane theory, in which the bending boundary layer 
is ignored such that the radial displacement is presumed to be constant valued. First, equation (B 1) 
gives a constant value for the axial stress resultant N x . Thus, the boundary condition given by 

equation (C2a) yields the result N x (x) = N L . Similarly, equation (5) yields T(x) = T 0 , where T, is 
an unknown constant. However, because this problem does not have a boundary condition 
specified for T (x), the problem is statically indeterminate. Next, the transverse-shear stress 
resultant Q x , that arises from variations in the bending stress resultant M x (see equation (B4)), is 

neglected in equation (B3) to get N e (x) = 0. Similarly, the derivatives of the radial displacement 

w vanish and, as a result, the rotation P° = 0 the bending strain k° = 0. The matrix constitutive 
equation given by equation (14) reduces to 
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\ W 

! R 


(C3) 
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In addition, equation (16a) reduces to 


A 22 j^ — A ■ i 


A, 


(C4) 


Substituting equation (C4) into (C3) yields 


( A 22 + Ai 2 ^ i 2 ) A 26 ^ 12 

/ U ')= A - 

* 11 *16 

A 12^62 ( A 22 + A 26 ^ 6 2 ) 

\ , ( n 22 
l V J 

*16 *66 


\T„/ 


(C5) 


Integration of these equations reveals that u(x) and v(x) are linear functions. The only possible 
way that the boundary conditions given by (Clb) and (C2b) can be satisfied is for v(x) = 0. Thus, 
equation (C5) can be expressed as (see equations (13)) 


/*u 

\ T „/ 


A„ A„l/ U '(*)\ . /a„\w 
A„A „]\0 / \A„/ R 


In addition, equation (C4) gives 


(C6) 


r 

U 



(Cl) 


Substituting equation (Cl) into the first of equations (C6) gives the result 


vy _ A 12 N l 

R A A -A 2 

-^ 11-^22 ^ 12 


(C8) 


Next, substituting equation (C8) into equation (C7), integrating, and applying boundary-condition 
equation (Cl a) gives the result 


u(x) 


-^22 N l X 
A 11-^22 ^ I 


Substituting equations (C8) and (C9) into the second of equations (C6) gives 


(C9) 


r-p _ A 22 A ]6 

Ml 


A n A 2 


A|2A26 n, 


(CIO) 


Approximations to the bending stress resultants are obtained by substituting v(x) = 0 and 
equations (C8) and (C9) into equations (17a), (82), and (88) to get 
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N, 


(C13) 


Because these bending stress resultants are constants, it follows that the transverse-shear stress 
resultants Q x and Q e are zero valued (see equations (B4) and (B5)). An approximation to N x0 is 

obtained by substituting T(x) = T 0 and equation (Cl 3) into equation (4). 
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Appendix D 

Approximate Solution for Uniform Pressure and Temperatures 

For this example, the boundary conditions are specified as 


u(0) = 0 

(Dla) 

v(0) = 0 

(Dlb) 

w(0) = 0 

(Die) 

|3°(0) = 0 

(Did) 

u(L) = 0 

(D2a) 

V (L) = 0 

(D2b) 

w(L) = 0 

(D2c) 

P° X (L) = 0 

(D2d) 


at x = L. The applied mechanical load is given by q„M = q° and the applied thermal load is given 
by ©o(x) = H 0 and 0,(x) = G 0 . For these loading conditions, equations (lb) and (7b) yield 
N,(x) = 0 and T,(x) = 0. Likewise, equations (100) reduce to 


Nl h (x) 

T T Th ^ Th 

= H 0 n x + G 0 c x 

(D3a) 

n:( X) 

— H o n 0 + G 0 c e 

(D3b) 

M™(x) 

= H 0 c x + G 0 m x 

(D3c) 

^Th 

T (X): 

_ Th _Th 


= H o n x0 + G o c x0 

(D3d) 


In addition, equation (B17) yields 
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NieV) 

x t Th Th 

= ^O n x0 ^O C x0 

(D3e) 

M?(x) 

x t Th Th 

= H 0 c e + G 0 m B 

(D3f) 

M™(x) 

= H 0 c xe + G 0 m xt) 

(D3g) 


An approximate solution is found by applying the equations of linear membrane theory, in which 
the bending boundary layer is ignored such that the radial displacement is presumed to be constant 

valued. First, equations (1) and (7) give constant values for the axial stress resultant , N x (x) = N 0 , 

and the torsional stress resultant, T(x) = T 0 , where N 0 and T 0 are unknown constants. Because this 

problem does not have any stress-resultant boundary conditions, the problem is statically 
indeterminate. Next, the transverse-shear stress resultant Q x , that arises from variations in the 
bending stress resultant M x (see equation (B4)), is neglected in equation (B3) to get 

N„(x) = q°R. Similarly, the derivatives of the radial displacement w vanish and, as a result, the 

rotation P° = 0 and the bending strain k° = 0. The matrix constitutive equation given by equation 
(14) reduces to 


/ u '( x )\ 
\ v '( x ) / 


*n *i6 /n 0 + N x (x)\ ^ j $ l2 \ w 

*, 6 *66j|f 0 +f Th ( X ) / \^62 / R 


In addition, equation (16a) reduces to 


(D4) 


. w 0_ Th , , 

22 — q n B + N e — A 12 U — A 26 V 


(D5) 


Substituting equation (D5) into (D4) yields 


( A 22 + A 12 ) A 26 ^ 12 


■^ 12^62 ( A 22 "*■ ^ 26 ^ 62 ) 


l/“'l 

A r 

ii 

> 

to 

*11 *16 

Jl v J 


*16 *66 


( N “ + 5 MHu; R+ Nr) 

lT 0 + T J l ^62 J 


(D6) 


Integration of these equations reveals that u(x) and v(x) are linear functions. The only possible 
way that the boundary conditions given by (Dla), (Dlb), (D2a), and (D2b) can be satisfied is for 
u(x) = 0 and v(x) = 0. Thus, the left-hand-side of equation (D4) vanishes and the resulting 
equations are expressed as (see equations (13)) 
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/No\ _ / A 12 \ w _ /n x \ 
\To/-\A 26 /R-\f T 7 


(D7) 


In addition, equation (D5) reduces to 


A 


W 

22 R 


= q°R + N 


Th 

0 


Substituting equation (D8) into (D7) gives the results 


/No\ 

\To/ 


(q:R+N:w Al2 \ ( N n 

l A 22 )\Aj-\f h l 


(D8) 


(D9) 


Approximations to the bending stress resultants are obtained by substituting u(x) = 0, v(x) = 0, 
and equations (D8) into equations (17a), (82), and (88) to get 

M, = |^( q :R + ND-M; h (DIO) 

-^22 

M,=|£( q :R + Nr)-M; com 


m,.=|^(c:r + n:)-m; <di2> 

22 

Because these bending stress resultants are constants, if follows that the transverse-shear stress 
resultants Q x and Q e are zero valued (see equations (B4) and (B5). An approximation to N x0 is 

obtained by substituting T(x) = T 0 and equation (D12) into equation (4). 
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Appendix E 

Approximate Solution for Linearly Varying Pressure and Temperatures 


For this example, the boundary conditions are specified as 


N x (0) = N 0 

(El a) 

_H> 

o 

II 

(Elb) 

o 

II 

o' 

(Ele) 

M x (0) = 0 

(Eld) 

at x = 0 and 

u(L) = 0 

(E2a) 

II 

o 

(E2b) 

w(L) = 0 

(E2c) 

M x (L) = 0 

(E2d) 

at x = L. The applied pressure load given by q„M = q! + q x and the applied thermal load is given 
by ©o(x) = H 0 + H, x and 0,(x) = G 0 + G, x. For these loading conditions, equations (lb) and 
(7b) yield N,(x) = 0 and T,(x) = 0. Likewise, equations (100) reduce to 

N x h (x) = (H 0 n x h + G 0 c x ") + (H,n x h + G.cDx 

(E3a) 

N e (x) = (H 0 n„ + G 0 c e ) + (H,n 9 + GjCe )x 

(E3b) 

M x (x) = (H 0 c x + G 0 m x ) + (H,c x + G;m x jx 

(E3c) 

T (x) (H o n x0 + G o c x0 J + + GjC x0 Jx 

(E3d) 


An approximate solution is found by applying the equations of linear membrane theory, in which 
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the bending boundary layer is ignored such that the radial displacement is presumed to vary like 
the pressure and temperature fields; that is, 


w(x) = W 0 + W, x (E4) 

First, equations (1) and (7) give constant values for the axial stress resultant , N x (x) = N 0 , and the 

torsional stress resultant, T(x) = T 0 , where N 0 and T 0 constants determined from applying 
boundary-condition equations (El a) and (Elb). Applying these boundary conditions yields 


N x (x) = N 0 = N 0 


T(x) = T 0 = T 0 


(E5a) 

(E5b) 


By using equations (E3), (E4), and (E5), the matrix constitutive equation given by equation (14) 
is reduced to 


\v'( 




, 

\ = 

*11 

*16 

/ 

*16 

*66 





/N„ + H„n™ + GdC™ + (H,n™ + G,c™)x\ / \ W, + W, x 

' ~ ' r (E6) 


\ * _Th _ Th ( Th _Th\ / ' \ - / 

\T 0 + H o n x0 + G o c x0 + + G ^^ Jx J 62 ) 

Integrating these equations and applying boundary-condition equations (E2a) and (E2b) gives 


/u( x )\ _ 

\v(x) / 


^11 & \ 


'(N 0 + H 0 nl h + G 0 cl h )(x - L) + + Q.c^fx 2 - L 2 )\ 

(f 0 + H 0 h x g + G 0 c x g)(x - L) + i(H,rC + G,c xe h )(x 2 - L“) / 

ft I2 \ 2W 0 (x - L) + W, (x 2 -L 2 ) 

+ \W 2R 


(E7) 


The remaining part of the solution is the determination of the constants W 0 and W r These 
constants are found by first noting that P 0 appearing equation (40), and given by equation (29b) 
reduces to the constant 


Po = 


[(A 12*11 + A 26 *i 6 )N 0 + (A | 2 * 1 6 + A 26 *66)To] 

RDu e 


(E8) 


Likewise, equation (102a) reduces to 


P,(x) = Y 0 + Y, x 


(E9) 
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for the given loading conditions, where 


Y n = 


l r ° , H o 


D,, c 


[^In ( n e n x (A 12 « u + A 26 A| 6 j n xe( A 12 *16 + A 26 « 66 j j 


+ R ^ Ce ^26*16) ^xe(A 12 5 16 + A 26 S 66 ) ) 


Y, — [q n + ti (n 0 n x (A 12 5 h + A 26 S I6 ) n x0 (A 12 5 16 + A 26 5 66 j) 


D,, e 


R 


+ (^0 Cx (A 1 2 A 1 1 + A 26 « 16 ) c x0 (A 12 * 16 + A 26 5 66 )) 


(ElOa) 


(ElOb) 


Substituting equations (E4), (E7), and (E9) into equation (40) yields 


W n = 


Pq + Yq 

4T 2 


(Ella) 



(Ella) 


Approximations to the stress resultants N 0 (x), M x (x), M 0 (x), and M x0 (x) are determined by 
substituting equations (E4) and (E7) into equation (16a), (17a), (82), and (88), respectively. 
Approximations to the transverse-shear stress resultants are then obtained by substituting the 
resulting expressions for M x (x) and M x0 (x) into equations (B4) and (B5), respectively. An 
approximate expression for N x0 is obtained by substituting equation (E5b) and the expression for 
M x0 (x) into equation (4). 
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Appendix F 

Exact Solution for Sinusoidal Pressure and Temperatures 

For this example, the boundary conditions are specified as 


u(0) = 0 

(Fla) 

v(0) = 0 

(Fib) 

w(0) = 0 

(Fie) 

M x (0)=0 

(Fid) 

N X (L) = 0 

(F2a) 

T(L) = 0 

(F2b) 

w(L) = 0 

(F2c) 

M x (L) = 0 

(F2d) 


at x = L. It should be noted that these boundary conditions are for simply supported ends rather 
than the clamped ends considered in the corresponding example presented herein. The applied 

mechanical load given by q n (x) = q“ sinf 3 ^ ) and the applied thermal load is given by 

0„(x) = H s sin(^) and ©,(x) = G s sin(^). For these loading conditions, equations (lb) and (7b) 

yield N,(x) = 0 and T,(x) = 0. Likewise, equations (100) reduce to 


Nl h (x) 

= (H s nl h + G s cl h )sin(^) 

(F3a) 

NlV) 

= (H s n"’ + G s c™)sin(^) 

(F3b) 

M x rh (x; 

) = (H s c x h + G s m x h )sin( 5 jj) 

(F3c) 

^Th 

T (x) 

= (H s m e h + G s c xe h )sin(^) 

(F3d) 
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In addition, equation (B17) yields 


z 

ii 

(H s n x r e h + G s c x0 h )sin(^j 

(F3e) 

m:(x) = 

(H s c™ + G s m™)sin(^) 

(F3f) 

M™(x) = 

(H s c™ + G s m™)sin(^) 

(F3g) 


Applying boundary-condition equations (F2a) and (F2b) to equations (la) and (7a) gives N 0 = 0 
and T 0 = 0. Thus, N x (x) = 0 and T(x) = 0. Next, the radial displacement is presumed to be given by 


w(x) = Wsin(^) 


(F4) 


where W is a constant to be determined. By using equations (F3) and (F4), equations (14) 
become 


\v'(> 




, 

u/ 

* 11 

*16 

/ \ 

*16 

*66 

V 




/H s nr + G s c:\ 

\ tt — Th ^ _ Th ( 

s^xB + G S C X0 J 


- If l f 12 ) + 5l2 ) \ s i n ( ) 

R \t* / L \* 62 // [L 1 


(F5) 


Integrating these two equations and applying boundary-condition equations (Fla) and (Fib) gives 


/ u ( x )\ , 
\ v W/ | 


*11 a l6 

*16 *66 


\H s n 


Th Th 

+ G s c 


/ R \M 


I HP? + G s cr\ w f *i 2 \ + W/ i* \ I LG _ / 

< Th >+ " ' " ' \^ 62 f jH \ 


JCX 


(F6) 


Substituting equations (F3)-(F5) into equation (17a) gives 


M x (x) = M x sin(^) 


(F7a) 


where 


— ff s ^B u (s n n x + « 16 n x6 ) + B 16 (s 16 n x + S 66 n x6 ) c x j 

+ GjBujSucr 1 + *i«Crt) + B l6 (5 16 c x h + 5 66 c™) - m x h ] 


+ 


W 

R 


Rjt~ 


D 1 1 + B 12 + B j , 


. Rjt . 

+ c. 


+ B, 


. RjF . 
+ ^2 


(F7b) 


Inspection of equation (F7) reveals that boundary-condition equations (Fid) and (F2d) are also 
satisfied. 
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The remaining part of the solution is the determination of the constant W. This constant is 
found by first noting that P 0 = 0 in equation (40) and equation (102a) reduces to 


P,(x) = Y sl sin(f^) (F8) 

for the given loading conditions, where Y sl is given by equation (102e). Substituting equations 
(F4) and (F8), and P 0 = 0, into equation (40) yields 


W = 


y s1 l 


ji 4 - 4S nV + 4TX 4 


(F9) 


The remaining stress resultants are determined by substituting equations (F3)-(F5), (97), and 
(B17) into equation (16a), (82), and (88). The results are: 


N e (x) = N e sin(^A) 


(F10a) 


with 


Ne — H s[a 12 («nn x + 5| 6 n x6 j + A 26 (5, 6 n x +5 66 n x0 j n 0 j 

■p G s [A 12 (* n c x + 5, 6 c x0 j + A 26 (a 16 c x + S 66 c xe ) — c e j 


W 

R 


Rnf 


B i2 + A 22 + Ai 2 | 12 


, Rjt . 

H — c 

L 


12 I + A 26 i c* 62 


. Rji . 

+ -j-r *6 


(FlOb) 


M 0 (x) = M e sin(^) 


(FI la) 


with 


M 0 — H s [B 12 (« u n x + S 16 n x0 ) + B 26 (s 16 n x + S 66 n x0 ) c 0 j 

+ G s [b i2 (5 u c x + « i6 c x0 ) + B 26 (« 16 c x +5 66 c x0 )-m 0 ] 


+ 


W 

R 


Rjt' 


Dp + B , 2 + Bp 


+ ^ I 


+ B , 


. Rjt . 

+ ~2 ~ if, 


(FI lb) 


M x0 (x) = M x0 sin(^) 


(FI 2a) 


with 


M x e - H s [B 16 (<t u n x + a 16 n x0 j + B 66 («i 16 n x + <*66^xe ) ^xe] 


+ 


G s [B 16 (s n c x + « 16 c x0 j + B 66 (s, 6 c x +5 66 c x0 j m x0 ] 


+ 


W 

R 


Rjt 


Di6 + B 26 + B | 6 | if |2 


. Rjt _ 

+ — r 


+ B, 


. Rnf - 

+ — r 


(FI 2b) 
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Likewise, equations (B4) and (B5) give 


qa: 

1 = ^ M x cos(f^) 

(FI 3) 

Qe( x ) 

= £ M xfl cos(^) 

(FI 4) 


An expression for N x0 is obtained by substituting T(x) = 0 and equations (F12b) into equation (4). 
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